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KEY   TO    SYMBOLS 


Y        V  =   Displacements 

a'      3 


"a'  \ 


>>,  y,  ^^ 


=  Displaceir.ent  components 

=  Cartesian  coordinates 

^'^  2,"  ~    Coordinates,  defined  bj  equation  (1) 

i         f  =  Thicknesses  of  core  and  faces,  respectively 
^c'  'f 

'      "  =  Superscripts  indicating  the  lov/er  and  upper 
face  quantities,  respectivelj' 

t  =  (t   +  t  )/t 

c    f    c 

y  -  Total  potential 

\\        \j  =  Sti-ain  ener[;ies 
c'   i 

^\'   V,'   V/  =  W'ork  performed  by  external  forces  and  moment: 
q '   N '   M 

c^  p;j  Y,  P-  =  Indices  taking  on  values  x  or  y 

i   j^  k  =  Indices  taking  on  values  x,  y  or  z 

e    e  ^.  e   ,  e, ,   =  Components  of  strain 

ij   .orp'   a.-i  -iA 

-f     T  ^   T  T     -  Components  of  stress 

ij    (up        gS  33 

G  =  Shear  m.odulus  of  t]ic  core 
c 

vi  =  Poisson'  s  ratio  of  the  face 
f 

E  -  Young' s  modulus  of  the  face 

p  =  Bending  rigidity  of  the  face 

A  . 

5  =  Kronecker  delta 


=  [(l-v)/2];5        6.      +6        0.      +[(2vJ/(l-vJ]5      6 


YM^ 


N 


^   Pre-huckliiig   axial    forces   per  unit    length 


q  -   Lateral  load  per  unit  area 

Mo  =  Externally  applied  moment  per  \ini  t  length 

-dp 

Q  -   Transverse  sheai-  forces  per  unit  length 

a 

l\   .  =  Bending  moments  per  unit  length 

N  r>  -  Axial  forces  per  unit  length  corresponding 

^  to  small  deformations  during  buckling 

(n)  =  Bracketed  index  indicating  physical  compo- 

nent of  a  tensor 

I 
f   ]  =  Chri  stof f el  symbol  of  the  second  kind 

'■nq-' 

X  =  General  coordinate  variable 

n, X,q , s , t ,m,p       =  Indices  taking  on  values  1,  2  and  3 

Xm  ,,  .  •  . 

e.g.  -  Metric  tensor 

^   '  ^Xm 

X  ,  Y  =  Arbitrary  tensors 

n   n 

r,  G  =  Polai-  cooi-dinates 

2 

V  =  Laplacian  operator 

u  =  Pre-buckling  lateral  displacement  of  m.iddle 

plane  of  faces 

u  =  Latei'al  displacement  of  middle  plane  of 

faces  for  axisymmetric  equations 

N. ,  N  =  Compressive  axial  forces  per  miit  length 

applied  to  inner  and  outer  edges,  respectively 

a,  b  =  Inner  and  oviter  radii  of  annular  panel, 

respectively 

P  =  a/b 

D,  E  =  Defined  by  equation  (65) 

dw 
^  ""  dr 

F,  G,  H,  I  =  Defined  by  equation  (72) 


A,  B  =  Defined  by  equation  (73) 

Tl  =  r/b 

N,  R,  Q             =  Defined  by  equation  (94) 

§  -  Tl  -  1 

A  =  Coefficients  of  power  series 

p,  =  Buckling  coefficient  for  uniform  axial 
stress  distribution 

C-  =  Arbitrary  constant  of  integration 

J  ,  Y  =  Bessel  functions  of  order  one 

D  =  Radius  of  convergence 
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The  buckling  of  an  annular  sandwich  panel  is  investigated 
using  the  theorem  of  minimum  potential  energy.   Governing  equations, 
dei'ived  in  cartesian  coordinates,  are  transformed  into  cylindrical 
coordinates  by  means  of  covariant  differentiation.   Considering  t]ie 
faces  to  be  membranes  and  assuming  an  axisyiiimetric  buckling  m.odc, 
the  equilibrium  equations  are  uncoupled  through  the  application  of 
an  improved  technique. 

For  a  clamped  outer  edge  and  "slidei-"  inner  edge  a  sej-ies 
solution  is  applied  to  tlie  general  problem  of  ?'adial]j'  varying 
in-plane  stresses.   Critical  stresses  ai'e  plotted  versus  ratio  of 
inner  and  oxiter  radii.   The  fifth  approximation  is  shovm  to  yield 
acceptable  results. 


CHAPTER  I 


INTRODUCTION 


A  sandwich  panel  is  defined  as  a  three-layei'  panel,  consisting 
of  tv.-o  thin  ontci-  layers  of  high-strength  material  between  wliich  a 
thick  layer  of  low  average  strength  and  density  is  sandwiched.   The 
two  thin  outer  layers  are  called  faces,  and  the  interniedi  ate  laj-er 
is  the  coi'e  of  the  panel  [1], 

Among  tlie  main  advantages  of  sandwich  construction  are: 
a  high  rigidity  to  weight  ratio,  good  thermal  and  accoustical  insula- 
tion, and  case  of  mass  production.   Some  examples  of  core  materials 
are  balsa  wood,  cellulose  acetate  and  synthetic  rubber.   However,  in 
more  recent  times  tliin  foils  in  the  foi-m.  of  hexagonal  cells  perpendic- 
ular to  the   faces  liave  been  employed.   Depending  upon  the  intended 
application,  faces  may  be  construcied  of  aluminum  alloys,  high-strength 
steel,  etc. 

The  practical  importance  of  sandwich  construction  came  into 
prominence  with  the  advent  of  the  airci'ait  and  space  industries. 
With  the  need  for  lighter,  stronger  and  more  stable  structural  com- 
ponents, great  emphasis  was  placed  upon  the  design  and  a.nalysis  of 
workable  sandwich  panels. 


Numbers  in  ^rackets  refer  to  the  Bibliography, 
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The  essential  difference  between  the  analysis  of  single-layer 
panels  and  that  of  sandwich  panels  is  tliat  the  shear  deformation 
associated  with  the  cere  of  a  sandwich  panel  may  not  be  neglected. 
Moreover,  initially  plane  sections  no  longer  remain  plane  during  bend- 
ing, and  the  existing  plate  theories  [2,3]  require  extensive  modifi- 
cations. 

Numerous  authors  have  contributed  to  the  development  of 
matliematical  theories  describing  the  behavior  of  rectangular  sandwich 
panels.   Two  of  the  more  noted  of  these  are  Hoff  [4]  and  Eringen  [5]  wlio, 
early  in  the  1950' s  illustrated  a  concise  and  straightforward  approach 
to  tlio  problem  using  the  theorem  of  minimum  potential  energy.   In  1960 
Qiang  and  Fbcioglu  [6]  inti-oduced  continuity  conditions  for  displace- 
ments across  the  interface  of  two  adjacent  layers.   This  modification 
has  been  shov.ri  [7]  to  contribute  appreciably  to  the  accuracy  of  the 
derived  equations,  wliile  introducing  no  additional  mathematical  compli- 
cations. 

In  comparison,  cii'cular  sandwi  cli  panels  liave  been  tlie  subject 
of  relatively  little  investigation.   In  1949  Eric  Reissner  [8],  neg- 
lecting the  bending  rigidity  of  tlie  faces,  solved  tlie  pj'oblemi  of  a 
cii'culai'  sandwich  panel  subjected  to  axisyimnetric  transverse  loading. 
Later  Zaid  [9]  included  tlie  effects  of  the  bending  rigidity  of  the 
face  layers.   Huang  and  Ebcioglu  [10],  using  a  teclmique  similar  to 
that  employed  by  Zaid,  recently  investigated  the  axisynmietric  buckling 
of  a  circular  sandwich  paiiol  subjected  to  uniform  axial  compression. 
In  t}ieir  final  results,  the  faces  were  treated  as  momljranes. 


Prioi-  lo  the  advent  of  sanchvich  construction,  the  stability  of 
circular  and  annular  single-layer  panels  were  extensively  investigated 
[3].   Jleissner  [11]  analyzed  tlie  axisyiimietric  buckling  of  a  single- 
layer  annular  panel  subjected  to  uniform  compj-ession  along  the  outei' 
boinidary.   In  this  vork,  the  inner  boundary  was  considered  free  and 
the  outer  boundary  simiDly  svippoi'ted  or  clamped.   Olsson  [12]  extended 
Meissncr's  analysis  by  considering  tlic  outer  boundary  clamped  and  the 
iniier  boundarj'  "slider."   Sucli  a  condition  could  be  approximated  by 
allowing  a  shaft  or  rigid  cj-linder  to  occupj'  tlie  central  hole  (see 
Figure  2) . 

Prompted  by  tlie  foi-egoing  sequence  of  investigations,  the  main 
objectives  of  tlie  present  analysis  are:   (i)  to  parallel  the  sequence 
of  analysis  present  in  the  literature  of  single-layer  panels  by  inves- 
tigating the  axisynunoti'ic  buckling  of  annular  sandwich  panels;  (ii)  to 
acliicve  a  more  satisf actorj'  foi'mulation  of  tlie  theory  througli  the  use 
of  continuity  conditions  [6];  (iii)  to  modify  the  uncoupling  procedure 
introduced  by  Zaid  [9]  and  later  employed  by  Huang  and  Ebcioglu  [10], 
tlicreby  signi  f  icaiitly  reducing  the  complexity  of  the  uncoupled  equilib- 
rium equations;  and  (iv)  to  apply  the  boundary  conditions  employed  by 
Olsson  [12]  lo  the  present  work. 

For  the  initial  derivation,  the  set  of  basic  assumptions 

employed  in  the  picsent  anal3'sis  are: 

(Al)   The  effect  of  transverse  normal  stress  in  each 
layci-  is  negligible. 

(A2)   Tlie  core  undergoes;  shear  deformation  only. 

(A3)   Displacements  in  each  layer  ai-e  linear  functions  of 
distance  from  Ihe   median  plane  of  the  layer. 


(A4)   The  median  plane  of  the  core  remains  neutral  during 
small  deflections  of  the  panel. 

(A5)   Each  layer  is  homogeneous  and  isotropic. 

(A6)   The  core  is  attached  to  the  faces  securely. 

(A7)   Hooke' s  law  is  valid  throughout. 

(AS)   Local  buckling  of  the  panel  does  not  occur. 
These  assumptions  are  similar  to  those  iised  by  Hoff  [4]  and  Chang  and 
Ebcioglu  [13]. 

Equilibi'ium  equations  and  boundary  ccndi  Lions  arc  derived, 
using  tlie  above  assumptions  and  the  th.corem  of  minimum  potential  energy. 
These  equations  are  then  transformed  into  polar  coordinates  thi-ough  the 
application  of  tensor  analysis.   In  order  to  solve  our  problem,  the 
resulting  equations  are  then  simplified  further  using  the  following 
add  i  1  i  on  al  as  su)iii:>t  i  ons : 

(Dl)   Bernou]  li-Navier  lij-potliesis  is  valid  foi'  the  faces. 

(B2)   Faces  are  considei-cd  to  be  membranes, 

(D3)   Annular  sand\^•ich  panel  subjected  to  uniform  axial 
compression  buckles  axisj-mmetrically. 

The  first  assumption  is  valid  wlien  the  sliear  deformation  of  the 

face  layers  is  negligible.   Kim  [7]  lias  sho\ni  this  to  be  true  for  most 

practical  applications.   The  second  assumption,  which  has  been  employed 

by  many  authors  [8,10],  is  discussed  at  length  by  Plantema  [1].   Finally 

the  assumi)tion  of  axisymmetric  bvickling  cannot  be  rigorously  verified 

without  solving  the  more  difficult  problem  of  unsyraraetric  buckling. 

However,  such  an  assumption  v.ould  seem  I'easonable  since  Olssoji  [12]  h.as 

proved  its  validity  for  the  analogous  probleji.  of  a  siogle-layer  panel 

witli   siifiilr.r  boundary  coi_di  tions. 


When  comparing  the  present  work  with  existing  and  subsequent 
theories,  the  words  panel  and  plate  may  be  considered  interchangeable. 
Numerical  computations  are  carried  out  bj'  a  desk  calculater. 


CHAPTER  II 
BASIC  CONCEPTS 

■^  •   r)cscrlption  of  Sandwich  Panel 
End  Notations 

An  element  of  tlie  sandwich  panel  to  be  considered  in  this  work 

is  sIiouTi  in  Figure  1.   For  convenience,  three  coordinate  systems  are 

defined  having  common  x-  and  y-coordinates ,  and  transverse  coordinates 

related  as  follows: 

t   +  t. 


(1) 

t   +  t^ 
//  _      _c f 

Z   -  Z  4  -    2 

wlierc  t   is  the  tliickncss  of  the  core  and  t„  is  tlie  thickness  of  the 
c  f 

face  layers  whicli  are  identical.   In  the  notational  system  employed 
here,  a  single  prime  denotes  quantities  related  to  the  lower  face, 
a  double  prime  those  related  to  the  upper  face,  and  vmprimed  quantities 
are  used  for  the  core  as  ivcll  as  for  general  relations.   Subscripts 
c  and  f  denote  quantities  related  to  the  core  and  face,  respectively. 
The  total  tliickness  of  the  panel  is  assumed  to  be  very  small 
in  comparison  with  the  lateral  dimensions.   Also,  each  layer  is  con- 
sidered to  be  isotropic  and  h'')mogencous,  witl:  the  pi'operties  of  the 
faces  characterized  by  E   and  V  ,  and  those  of  the  core  characterized 
by  G^. 


3     z 


3     i' 


X       - 


-X       -- 


Fic-ure   1.      Element   of    Sandwich   Panel 


Arbitrary  transverse  load  q,  axial  forces  N  ^,  and  external 

— O'P 

moments  H  g  are  the  only  externally  applied  forces  or  moments  which 
may  contribute  work  to  tlie  system. 

For  the  initial  dei-ivation  and  formulations,  indicial  notation 
will  be  employed.   Greek  indices  0! ,    p,  y   and  |j,  will  take  on  values  1  or 
2,  while  Latin  indices,  wlien  employed,  will  take  on  values  1,  2  or  3. 
Any  repented  index  denotes  summation,  wliile  any  number  of  indices 
preceded  by  a  comma  indicates  partial  differentiation  witli  respect  to 
the  coordinate  variable  represented  by  those  indices. 

2.   Displacements 

In   the  following  relations,  w]iich  are  valid  for  small  displace- 
ments, it  will  be  assumed  that  within  each  layer  of  the  sandwich  panel 
plane  sections  are  preserved,  although  not  necessarily  perpendicular  to 
the  deflected  middle  surface.   This  will  constitute  an  approximation 
since  the  presence  of  transverse  shear  suggests  a  non-linear  variation 
of  the  longitudinal  displacements  through  the  thickness.   Also,  while 
the  panel  is  geometrically  syimiietric  with  respect  to  the  middle  surface, 
the  displacement  relations  will  not,  in  general,  reflect  this  syiimietry, 
since  unequal  bending  moments  may  be  externally  applied  to  each  face 
layej'.   llieref  ore,  wi  tli  the  aid  of  equation  (1),  our  displacements  can 
be  v.ritten  in  the  following  form: 

V^(x,y,z')  =  u^(x,y)  4  z'.;^'^(x,y)  V^(x,y,z')  =  w(x,y) 

VQ,(x,y,z)  r.  zv^,^(x,y)  y^(.^,y,y)    =   w(x,y)       (2) 

V^(x,y,z")  =  u^(x,y)  +  z/^.^(x,y)  V^(x,y,z")  :.  w(x,y) 


where  (V  ,V  )  is  the  displacement  of  a  generic  point  in  the  sandwich 

panel;  u   represents  the  displacement  in  the  xy-plane  of  a  point  lying 

in  the  median  plane  of  eitlier  face;  i},  represents  the  angle  that  the 

normal  to  the  median  plane  of  each  layer  rotates  v.hen  the  planes  are 

deflected;  and  w  is  the  transverse  deflection  which  is  assumed  to  be 

constant  througli  the  entire  thickness. 

In  Older  tf  ensure  continuity  at  the  interface  of  any  two 

adjacent  layers,  the  following  conditions  must  be  imposed  on  equation  (2) 

t.      t 
/  _  _f  I  /  _  __c  I 

(3) 

"0^  "^  2  ^'a  ~  "    2  ^u 

Using  (3)  to  eliminate  ■!■      and  u   from  (2),  our  displacement 
relations  become: 


\ 


\'      ^  u,  +  z  V 


(4) 


V^  -  -U^  +  z  A   +  (A   -  4  ) 

a     ct  'Or    2  ^cy   *» 


V3.  w 


3.   Strain-Pi s]>l  acement  Relations 

Linear  strain-displacement  relations  are  defined  by  [11]: 

®C^3  "  2-^cy,?  ^  ^'g,»^ 

^C.3  =  l^^'c.,3  ^   ^.e.)  <^> 


e   --  V 
33     3.3 
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Using    (4)    in   conjunction   with    (5),    the    sti—iin-displacement 
relations   for    each    layer   becomej 

e'      =    0 
33 


-   u^         4    — (y„         -    Vr^       )    +    z    w    ^       ] 


<3   =   Ik  ^   '\a^ 


e        -    0 
33 


4.       Stress-Strain   Relations 


(6) 


Tiie  generalized  Hooke' s  law  for  a  liomogeneous  isotropic  body 
can  be  v.:-itten  in  tlie  following  foi^n  [15]: 

T.  .  --.  -^[e.  .  +  -^  6.  .  e,  ,  ]  (7) 

where  the  Latin  indices  i,  j  and  k  take  on  the  values  1,  2  and  3. 
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\Vlien  the  transverse  normal  stress,  T   ,  is  neglected,  (7)  may 


be  rewritten  as  [15]: 


©•3         1-V      o;^         1-v      £<-3      YY 


T  =:     0 

33 


or, 


^       A 


^p      i-v^    ^f^Y^^   Y^^ 


(9) 


whei'e, 


A  „        =  2:::!:i(6        6„     +  6        6,     +  ;^  5    ,.   6      ) 


wi  til , 


A,         =A„  =A_         =A        „ 

C.cy[i  ^oy[i  a?i[Ly  y[uyp 

If  the  lateral,  in-plane  stresses,  T  „,  are  also  neglected, 
(9)  reduces  to: 

T   =-^e^=2Ge^  (10) 

03    1+v   0-3       Q'3 

Assumption  (Al)  implies  that  equations  (9)  may  be  used  for 
the  face  layers,  while  assumption  (A2)  suggests  the  use  of  equation  (10) 
for  the  core. 
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5.   Stress  Resultants 

In  order  to  reduce  the  tliree-dimcnsional  elasticity  problem  to 
a  two-dimensional  one,  the  following  stress  resultants  are  defined: 


^3  ^^ 


*f 

2 

^f 
2 

dz '     ; 

><;? 

s 

*f 

2 
2 

\ 

*f 

z'dz' 

^  ";3 

_^f 

2 

^3  ^  ^"  '  ^i^> 


t  t  t 

f  f  c 

2  2  2 

Q  s     T  ^  ciz   ;   Q   F     t\  dz'  ;  Q  E      T  ,  dz 

t  t  t 
i_                                       f               c 

~  2  "2  2 


where  M,  N,  and  Q  denote  internal  moments,  in-plane  stress  resultants, 
and  shear  resultants,  respectively. 

Expressions  for  stress  resultants  in  terms  of  displacements 
are  obtained  hy   substituting  strain-displacement  relations  (6)  into 
equations  (9)  and  (10),  yielding  relations  for  the  faces  and  core, 
respectively.   Wlien  these  relations  are  then  substituted  into  equa- 
tions (11),  and  tlie  indicated  integration  performed,  the  following 
expressions  for  stress  resultants  and  moments  are  obtained: 


■^W  =  °f  ^«3Y^l  ^.^x 


m''   =  D   A  „    ii" 
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N  _  = A_ 


C^3    ^_^2  *  C'^YM-  Y'^^ 


Op 


f    f  ^  f  1 

=    -^  A  ^    i  -u'    +    -~    (t'         -  f '   ) 


(12) 


t  E 


//     f  f       // 
^a        2(l-iV',)  "^^(7     .o' 


Q  r.  2G   [u'  +  -(t  w   -  t.  ly')] 
O'     c   a       2  c  ,»    f  ^cy 


whore  the  bending  rigidity  of  the  faces,  D  ,  is  defined: 


^f  ^ 


f  f 


12(l-v^) 


chapt]:r  III 

DERIVATION  OF  EQUATIONS  IN  CARTESIAN  COORDINATES 

1 .   Total  Potential 

The   total    potential    consists   of    the   strain   energy    stored    in 
tlie  panel    dui'ing    a   small    def oinnation   and   the   total    work  pei'foi'med   on 
the   panel    by    the   external    forces    and   moments    dui-ing    that    deformation. 

In  general,    the    strain    energy   stored    in   the  body    is    defined 
[14]: 


U  r_-   i-   rrr  t         e        dv  (13) 

2  J^c>       ij       ij 


In  pa^'ticular,  as  a  result  of  assumptions  (Al)  and  (A2)  ,  the 
strain  enerq,ies  associated  witli  the  faces  and  core  are: 


<  -  I   rjr  C^'.    ^'r   +  2t'   e'j  dv 


U  =  i-  fff  (2T  „  e  „)  dv  (14) 


""l-Ull    '\^    V-^'^^3<3^  ^^" 


where  the  integration  is  performed  over  the  volumes  of  the  individ- 
ual layers. 
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Also,  since  the  only  externally  applied  forces  or  moments 
which  are  allowed  to  do  work  are  the  transverse  loading  function,  q, 
the  axial  forces  N  ^,  and  the  external  moments  M  ^,  the  work  performed 
on  the  panel  is  defined  bj': 

W  =  rp  qw  da 
a 

1  r-r  X. 


K   =  -  T  n  N  ,(V'  o   +  V'    +  w    w   )  da' 
a 

K^  "  -  7  rr  N  n(v"     +  V"    +  w    w   )  da" 


(15) 


W  =  r  (m'„  ('   +  M-''   (/')  n_ds 
s    ' 

where  the  integration  is  perfoi'med  over  the  ai-eas  and  boundaries  of 
the  individual  layei's. 

In  the  second  and  third  integral  expi'cssions  of  (15)  ,  w]iich 
are  analogous  to  tliose  employed  by  Timoshenko  [3]  and  Eringen  [5]  ,  it 
is  assumed  that  tlie  total  axial  load  is  evenly  distributed  between  tlie 
uppei'  and  lower  faces. 

Therefore,  the  total  potential,  whicli  is  defined  as 

1r  =  V'    +  V      +  U"  -  W  -  W'  -  \/'  -  W  (16) 

f     c     f     q     N     N     M 

may  be  rewritten,  using  (14)  and  (15),  in  the  following  form: 
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y  =   i-   fff    (t'      e'„    +    2t'  e'    )    dv' 

V 

+   Ini    (2T    .    e    ^)    dv  +   i    nj    (t\    e%    +    2t''      e\)    dv" 
^  V 


(m'      i''    +   M^'      i;,'")    n      ds  (17) 

-a3     cy       -a3    V       3 


2.   Theorem  of  Minimum  Potential  Energy 

The  theorem  of  minimum  potential  energj'  states  that  of  all 
displacements  satisfying  the  given  boundary  conditions,  those  whicli 
satisfy  t]ie  equilibrium  equations  make  the  total  potential  energy  an 
absolute  minimui;i  [14],   Tliercfore,  equilibrium  equations  and  boundai-y 
conditions  for  tlie  sandwich  panel  are  given  by  the  variational 
equation: 

6  y  =  0  (18) 

However,  before  applying  this  extremum  principle  to  equation 
(17) ,  we  note  that: 


'  <V3  <s'  -  '  I  - 


r   ■^■r        /    ~| 


f 

2E, 


— -  A  „    e  „  6  e 


1_V^   apyl-''   CrP     YH- 


2."  „  6  e^„  (19) 
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and   similarly, 


6    (t'      e'^)    =    2t'      6e',^  (20) 


6    (T    „    e   ^)    =    2T    _    6e   „  (21) 

Q'3      cy3    ,  Q'3        Q'3 


6    (T%    e"  J    =    2T%    6e%  (22) 

q3      qo  q'3         q'3 


6    (t"       e". )    .=    2t'^    6e''  (23) 

op       QP  Qp  Q^ 


Therefoi'e,    using    equations    (19-23)    in   conj\n:ic.tio3i   with 
displacement    relations    (4)    and    strain-displacement   i-elations    (6), 
the  variation   of    equation    (17)    yields: 

6y  =    IT    ^(ou      r    +    7.   6v      ,)    +    T    ,  (6((      +    6\v      )      dv 

V 

.u;.|        c'j  L  t.  2         .Q-       t         Q'  J 

V  c  c 

+         Ht   „      -   6u  +   z   6i^  +  (6^,  _   6-^        ) 

••'<';'    I.    GP    L  <^,P  CV.P  2  ^Q',3  ^Q,P    J    • 

v' 

+   t''        6  a'    +   6w  i-  dv" 

Qf3  L      »  jOj  -• 

-  ff  q   6\v  da  +   -    ff  N   .  r2u-        6\v  ^    +   ~   (6i|i'    „    -   6a"    ^)      da 

-  r        m'      6iV'    +   m"      6il'''      n-,    ds   =   0  (24) 
^      L-G'P      ^Q-        -cvS      ^o-J     p 

In   equation    (24)    the   syiimietrj'   of    the   stress    tensors    and 
strain   tensors   has   been   employed.      Also,    the   displaceme:its 


associated  with  N  „  arc  those  occuriing  in  the  median  planes  of  the 
face  laj-ei's. 

Using  relations  (11)  to  reduce  the  volume  integrals  In 
equation  (24),  and  applying  the  two-dimensional  divergence  theorer,; 
of  Gauss  to  the  resulting  area  integrals,  the  final  integral  equation 
is  obtained. 

As  an  example  of  the  above  procedure,  consider  the  first  Xai-ia 
in  the  integrand  of  equation  (24): 


t^/2 


da 


TJ/t'   6u'  ,  dv'  =  rr  J     t'   6u'  ,  dz' 

/  a  -t^/2        ' 

V  f 

=  ff  n'   6u'    da 
a 

.  n[(^\   6u')  ,  -  n'   .  6u'"]da 

=  f  n'  6u'    n   ds  -  rr  n'   „  6u'  da  (25) 

s  a 

After  applying  a  similar  procedure  to  the  remaining  terms, 
and  collecting  coefficients  of  similar  virtual  displacements,  a 
complete  description  of  tlie  sandwich  panel  in  the  equilibrium  state 
is  obtained  in  integral  form: 


6r 


rr  {  T-  q  -  Q  -  N  ^  w   _  -  q'    -  q"   "1  6w 

''^     <  L.                Of, a  -ap      ,cyp  a, a         o:,<yJ 

p2Q        ^  -  r    ■''f 

L  t     CiP  ,3  ap ,  pj  a  L  t  a         &f3 , 3 
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+    Q'    -   -^   n''   ^1  5ik'  +  r~  n"  ^  „  -  ^f'   „  +  q"  I  6f'|  da 
a         2   0-3,  f J  ^a       L2    cy^,3    oP.P    oj  ^d 

+  r  i  Q  +  N  „  w  .  +  q'  +  q"   6w 
^       I  \_  a       -op   ,p    a    cyj 

L  cc?>  o.SJ   3   \_6    -a^  q3    2   q-3    -la-gj   ^  P 

+   -  ~  N  ^  -  -^  n''   +  m''   -  k'   I  6-^''}  n   ds  ^  0  (26) 


3.   Equilibrium  Equations  and 
Boundary  Conditions 

If  v.e  apply  tlie  fundamental  lomma  of  the  calculus  of  varia- 
tions [15]  to  equation  (26),  we  obtain  the  following  equilibrium 
equations  in  terms  of  stress  resultants  and  moments: 


N  „  w   5  -  (Q   -I  0  +  Q  )    =  q 


2Q 

T^-<3.3  "^''aP.B^  ° 


t       ^        t 

t  ^   Q'     Q-3  ,  3   Q'     2    Q§  ,  3 


-^  n\  „  -  m''   .+  q"  =    0  (27) 

2   Qi::> ,  3     QP ,  P    <y 


and  boundary  conditions, 

■   r  I  Q  +  q'  +  q"  +  N  „  w  ,  I  6w  n   ds  =  0 

'J   L  a     cy  Q-  -Q'3   ,3J      Q' 

r   n',  -  n''  6u'  n^  ds  =  0 
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r-t  t      -, 

If  v;e  substitute  the  stress  resultant-displacement  relations 
(12)  into  equations  (27)  and  (2S) ,  and  collect  similar  terms,  we  obtain 
the  following  equilibrium  equations  in  terms  of  displacements: 

1                                S^f 
-  2G   [u'    +  -  (^  \^'     -  tj'       )]  -  N  „  w  a  -  -— i [2w 

c   a,cr'       2       c    ,(ya         fa, a  -a3  ,op        2(l4V  )    ,aa 

+    (^•'   +'i"   )]  =  q 
"a, a   'a, a 


4G     ^  t  E 

-— ^  [u'  +  -  (t  w    -  t^il')]  +  —   ■  ■■  A  ^    [-  4u'  ^   +  t  (*'  ^ 
t^        a        2   c  ,a  i^a''        2(l-v^)  ""'^"^'^  V.^^^    f  \,P^l 

-  f'    )]  =  0 


-2t,G  ,  t^E^ 

^C^/  1,,            j,'sn  ff     ri'          T 

— [u   +  -  (t  w    -  t    i;    )]  +          [((   +  w   ] 

t      a  2   c  ,Q'    IV  2(l4V^)  ^a          ,a 

C  1 

+  ^  A  -    u'  -  D.  A  .    [4^   „   -  3il'"  f.  ]  =  0 

2(l-v  )   ^'^"^^  '^'^'■^  '^'''^^    "i'^V-           Y.PM' 


-t^E  t  E 

—  \      A  „    u'  „   +  — i- i—-  [(/'  +  w   ]  -  D^  A  ,,    [-  3*'  . 

2(1-V^)   ^">'^''  '^'^^'    3'^^+^-^f>   ^    '^     ^  ^^v^^     *7.P^i 

+  4f'    ]  =  0    (29) 
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and  boundary  conditions, 


f    f 


f     ^2G    [u'    +  i   (t   w        -   t^^')]    +  N  p    w  .    +    ^,:    ^    ,    [2w 
0       Leo-        2        c    ,a  f^Q-  — ^        -•        "/-...    \ 


-QfP       ,P  2(l  +  v^)  ,« 


■+   a'    +   ilr"]  r   6w  n     ds  =      0 


/     \ ^A  [4u'         -   t    (^'  -   f      )]ku'    n„    ds   =    0 


i— -  ^    „    + — -  A„  —    (i)  -\k         )-tu 


+    D^   A   -.         \|f'         -   M   -,^    6<j('    n,,.    ds   =    0 


r 


t.  t„E„  r-t^ 


f   -^  N   .    -   — Xf      A   .        T-^-   (/         -    f      )    ..   t,   u'      ] 


s  2(l-v^) 


+    t)„    A   ^         i<  -   M   . 


lU      n,,    ds   rr    0  (30) 


Boundai-y  conditions  (30)  state  that  either  the  quantities 
contained  v.ithin  tlie  large  brackets  must  vanish  along  the  boundary, 
or  the  variation  of  the  displacement  r.iust  vanish  along  th.e  boundary. 

Compared  with  equations  (27)  and  (28),  equations  (29)  and 
(30)  constitute  a  formulation  tlirough  which  the  buckling  problem  can 
be  more  I'oadily  solved.   This  becomes  inmiediai  ely  evident  by  consider- 
ing the  classical  concept  of  instability  which  is  associated  v.'ith 
displacement  fields  rather  tlian  stress  fields. 
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4.   Simplification  of  Equations 

In  ordei-  to  simplify  equations  (29)  and  (30),  we  exclude  the 
possibility  of  externally  applied  edge  moments  by  specifying: 

m'„  =  m"   =  0  (31) 

-ap       —dp 

which  reduces  the  number  of  dependent  variables  from  seven  to  five 
since  (31)  suggests: 

il'    =  f  (32) 

a  a 

Also,  if  we  assume  that  plane  sections  remain  poi-pendicular 
to  the  deflected  middle  sui-facc  of  each  respective  face  (asssumption 
Bl) ,  we  have, 


(r   =  -  w  (33) 

a  .a 


whi  c]i  further  reduces  the  number  of  dependent  variables  to  three. 

A  final  assinnption  (B2)  which  simplifies  the  problem  by 
reducing  the  order  of  tlic  differential  equations  is  [1]: 


t  »  t.  (34) 

c     f 


whic]:  justifies  neglecting  the  bending  rigidity  of  the  faces,  D  , 
compared  witli  tlie  bending  rigidity  of  the  panel  as  a  whole. 

Conseqviently .  when  equations  (29)  and  (30)  arc  simplified 
throxigh  tjie  application  of  (31-34),  we  obtain  the  following  fc..;ui- 
lib.tium  equations: 
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2 

-t    t      G      (-—  u  +    £    V.'        )    -   N   Q    w  =    q 

c     c     t       cyjCv  .Q'Q'         -ap      ,ap 

c 

2t    E 
2G      [-^  u      +    £   w      ] ~  A  ^         u      ^      =    0  (35) 


and  boundary    conditions, 


•{2t    t      G        -^  +   -  V,-  +   N  o    w   r.  I    'Sv.-  n      ds   =   0 

"J        L  c      c  Lt  2      jqJ        -op       ,PJ  ex 

s  c  , 

r2t    E  . 

r     1 — ■— ^  A   .        u         ^    6u      n„    ds   ^   0  (36) 

s       I      .^ 

Equations    (35)    rn;d    (36)    can  be   rewritten    in   extended   form   as: 


-    i:    t      G 
c      c 


r  2  ^        ~\     ^  r  2  f        '] 

—  u  +    i    w  -    t    t      G        -—  11  -I    t    w 

Lt         X,:;  ,xxj  c      c  Lt„      y,y  ,yyj 


-   N        v.-  -    2X        w  -   N        w  =    q 

-XX      ,xx  -xy       ,xy        -^-y       ,yy 


T'  2  ^           "I        ^f^f 

2G     I  T-  u      +    t    w [d+V    )(u  +    u  ) 

c  It  X                 ,x  I        ,2                f         x.xx           y.xy 

"-    c  '    -        1-V                                '                 ^  ' 


+    (1-V^)(u  +   u  )]    -   0 

f        x,xx  x,yy 


;     r  ~  u      -1    f    -.v       1  -   —-4    [(l-v^)(u  -;-   u  ) 

c  Lt^,      y  ..yj        ,_^^2  i        y,yy  x,xy 


+     (1-V     )  (U  +1.1 

f       y,xx         y,yy)J   ^0       (o/) 
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and, 


+  N        w      r6'An      dv   =   0 
X        -xy       ,yj 


f    -^  2£     t     G  ~    4      -    W  +    N  V/ 

•J     L  c    c  Lt^         2       ,xj        -XX      , 

X   ,  pU  <^  — j  .^ 

f  J2f    1   G        -^  4    -  w  +   N        w        +   N        w      f   6w  n 

•J     L  c   c  Lt^        2      ,yj        -yy      ,y        -xy      ,xJ 

^     l~T-L^,x-^^f    "y,yj|^\"x^^=° 
f 

i    llT^   f"x,y   ^   "y,xV   ^^   \   ^>'   =   ° 

^x    2t^E_^ 
•^     ^-1    ,2         y,y  f       x,x  J       y      y 

f  ir—-—   [u  +11        ][    6u     n      dx  :^   0 

^     LlHV^         x,y  y,x  J  X      y 


dx  =   0 

y 


(38) 


1  <;  C  f 


and   primes    and   double   primes   have  been   oiiutted   without    confusion   as 
a   2"e5ui  t    of    equation    (31). 

Equations    (37),    which    ai-e   a   special    case   of    equations    (29), 
are   identical    with    those   derived   by    CTiang    and   Ebcioglu    [13]    if    thermo- 
elastic    effects    are  neglected. 

It    should   be  noted    tliat   when   simplification    (33)    is    introduced 
into   equations    (29)    and    (30),    equations    (,'3r->)    and    (36)    are  riot    imme- 
diately  pi'od^iced.       The   lasL    two   of    equations    (29)    niust    be   i-eturncd    to 
the   area   integral    of    equation    (26)    where,    as    a   result    of    (33),    tJicy 
become   coefficients   of    -5w      .       Conseouently ,    a    transformation   proce- 
dirre    similar    to    that    i  J  lu;.ii-at  en    in   eqi.a'.iou    (25)    yields    additional 
terms   v.hich    contribute   to    equatiois    (35)    and    (3G). 


CHAPTER  IV 


EQUILIBRIUM  EQUATIONS  AND  BOUNDARY  CONDITIONS 
IN  CYLINDRICAL  COORDINATES 


1.   Application  of  Covai-iant  Derivatives 

Let  us  consider  a  covariant  derivative  of  an  arbitrarj' 

covari ant  vector,  X  ,  [16]: 
n 

Sx 

X  I   =  — i'  -  {  ^}  X,  (39) 

n|q    ^  q    'nq-'  i 

.       In  (39)  y   is  a  general  coordinate  variable,  f   1  is  the 

nq-" 

Chrjstoffel  symbol  of  tlie  second  kind,  and  all  indices  take  on  values 
1,  2,  or  3.   Eui'tliei'niore,  when  a  vertical  line  precedes  any  nnniber  of 
indices,  it  indicates  covariant  differentiation  with  I'espect  to  the 
cooi'diiiate  variables  represented  by  those  indices. 

In  considering  the  second  covariant  dei'ivative  of  an  arbitrary 
vector  we  must  recall  tliat  in  order  to  sum  two  indices  we  must  have 
one  covariant  and  one  contravariant  index.   Therefore,  inti'oducing  the 
metiic  tensor,  g   ,  (39)  implies: 

X  ,^=  g-^^^X  ,   =  g'^f-^-  {  ■}  Xl  (40) 

n|     *=    njq       \_      q.        "-nq^  ij 

And  since  tlie  covariant  derivative  of  an  arbitrary  second  order 
mixed  tensoi',  Y  ,  is  [1  61: 
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26 


ox 


we  have,    fi'om    (40)    and    (41): 

'Sx"  ■-  Bx 

8X 

Sx 
Bx 


n|pt   -   %s      n|t   -   Ss   \^t  L^        ^^q 


r  m  ,   r    sq      ■  m  ^Q    r   ^^  -i    ^-  ~1  1  ,  .^^ 


In    cylinrlri  cal    coordinates,    the    Christoffel    symbols    and 
iTietric    tensoi'   take   tlie   following   form    [16]: 


^2-^         ^21^  1    '  '■22-'  ^ 


all    other    [       '(    =    0 

11  33 

^11    ~    ^33   "    ^        "    ^-        ~ 

1     *>  99  1       9 

529   =    <X    )"    ;  g        -     (—)''  (43) 

X 


Equations    (42)    and    (43)    imply: 


X    I     ,    -    X    I  (44) 

n  pt  n  tp 


Therefore,  -.vith  the  aid  of  (43)  and  (44),  (39)  and  (42)  m?y  be 
expanded  to  ji.eld: 
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SX 


11 


ax 


'  OX 


^1 1  2   "    ~2'  ~   "T  ^2 

'  dx  X 


^^'2    1 

^2|l    -   71    ~   "T\ 
'  Sx  X 


3    2 


Sx, 


Sx' 


3    1 


(45) 


and, 


S"x, 


1  11       ^,  1,2 
S(x  ) 


_^i    a!!2    1^   ^. 

'1J22  2,2  1,2^^      .    1    ~   -^1 

'  s(x  )        X    ax  Sx 


1    12 


2    11  .,    12  .1^2  1    -,    1 

'  d(x   )  (X   )  X      Sx 


'"''^  ,   1    ^\  1    ^\        ^^ 

^2l92   = 2~¥  '    ^^      ~2   +   >^      ~T  "    2X 

'  S(X   )  3x^  SX^  ^ 


a^x, 

1  ^''i 

1^2 

X     Sx 

+ 

^\ 

1     ^^2 

ax'sx" 

(X^^ 

1    .    1 

X     Sx 

'^\ 

2X^ 

2     S^2 

X  t\  ^  ^^2  1    ^1         ^ 

2112   ""         ].^    2  1    .     2   "•    ^      ^    1  1 

'  Sx   cX  X      Sx  Sx 
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d\ 

^3|ll    = 

3(X  ) 

^3|l2   = 

1     ^^3 
.1    -^.,2 

3X  3x    X  3x' 

^'^3      1  ^\ 
'      3<X  )        3x 

In  order  to  make  each  equation  of  (-15)  and  (46)  dimensionally 
homogeneous,  we  replace  the  right-hand  members  by  their  physical  com- 
ponents, X    ,  tlirougli  t}ie  relationship  [16]: 

X ,  .  =  -v'  g   X       (No  summation)  (47) 

(n)  n 

v.'here  a  bracketed  index  indicates  physical  component. 

Furthermore,  to  ensure  dimensional  compatibility  among  equa- 
tions, the  following  relations  are  iised  to  introduce  physical  compo- 
nents to  the  left-hand  members  of  equations  (45)  and  (46): 


\n)|(X)-^^^   -'  ^        \U 

I    nn  /~ li     I    pp"         ,  ,,„, 

^(n)I(X)(p)  ^  -^e   yg   Vg   X^i^^^   (Nc  summation)     (48) 

TJie  left-hand  members  of  tlie  resulting  equations  may  be 

interpreted  as  cartesian  comjionents.   Therefore,  after  replacing 

1         2 

X   by  r,  X   by  9,  X    by  u^ ,  X^^   by  u^,  and  X^^^  by  w  in  the  right- 
hand  members,  and  X    by  u  ,  X    by  u  ,  and  X    by  w  in  the  left- 
hand  members,  transfor;.;at  j  on  equations  relating  cartesian  and  polar 
coordinates  take  the  following  farm: 
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u     =   u 

y        e 


1  1 

y   y        r      6,9        r      r 


1  1^ 

xjy  "    r      r,e        r 


W;  =    W 


y       r    ,e 


X   XX  1- ,  rr 


'xjyy  ^  7  ^.ee  -  ~2  "e,e  "^  ?  ^,r  "  ~2  "r 


^Ixy   =    I  "r,re    "   "2   \,e    "   ^  ^G.r   "^    \  "o 


yjxx  G,rr 


V|yy  "■  3^"  ''e,6e  "^  "I  ^,e  '  F  "e.r  "  "T  "e 


V|xy    -    i   "e,r9    "   -2   "0,6   ^    ^  ^^ ,  r    "   "^  ^ 


.1  .■       _  2 


yy 


1     ^        ,  i  ,, 

'2   ^'',00   '    r    '\r 


(49) 
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In  equation  (-19),  r  and  9  are  polar  coordinates  defined  in  the 

usual  manTier,  while  u  ,  vi   and  \v  represent  displacements  in  tlie  direc- 
r    9 

tions  associated  with  tliis  new  coordinate  system. 

2.   Transf  ormat  j  on  of  Eqviations 

If  we  interpret  the  derivatives  appearing;  in  equations  (37) 
and  (3S)  as  covariant  derivatives  and  thereby  transform  them  into 
polar  coordinates  using  relations  (49),  we  arrive  at  the  following 
equilibrium  equations: 

-t  to  {AriA(ru).lu^  1.?.  V^w} 
c  c  Ll   1  r  or     r    r   6,9 J         J 
c 

-  ^r  ^rr  "  ^^-G  ^^   ".rO  "  ^2   ^e^  "  ^96  ^^  ".99  ^   I  \r^    ^    ^ 

r  r 

^^^f^V  r  2  1  /''f"^ 


J   — -  U   4  t  w — -   7   u   -  — T-  u 

cLt    r       ,1'J   -,   2  L    1-     2 
c  1 -V .  r 

r"-  ("9,r9  -^^,99^]=  ° 

r  2     t   1  'f ""f  r 

L.l^   t3    r   ,9j   i..,,2  L    f 

1  ^''~"'f'     1 

-  u    .)  +  -;i-  U   .,   --  0 

r   r,rG       ^  r,o_J 


(1+v^) 

-I 


—    t  E 
I     If  I  ?        1 

2GJ  ~  u„  4  £-  .,  J  _  -i-i  I  (l-v^)(v^  Ug  -  -'^-  Uq) 


1          1          ^^'-~^f^ 
+  (l-i-v^)(--  Ug  p^^  +  -  u..  ..J  + u__  .,  I  --  0     (50) 

r    ' 


and  boumlaj-y  conditions, 
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+  N  ^  —  w  ^h  6w  n   dG  ---^   0 
,  r   — r9  r 


y,     p-u    A   -, 

f  -^2tt  G   -^  +  -  w    +  N    w 
>>     l        celt     2   ,rj   — rr   , : 
"-  c 

6^2t  E  p  1     ~1^ 

r  -j —  j  u    +  V   (-  u    -h  -  u  )  \\6u      n   de  •-.  0 

•J  L   2  L  1' ,  1'    f   1   0 , 9    r   r  JJ   r   r 

f 

r  {r^  [-   u     -  --  u   +  u    I  }  6u   n   dG  --  0  (51) 

■   U+v  Lr   r,G    r   G     6|rJ  J    0   r 


2  2 

2        S  IB    13 

where     V  =  ---  +  -  ---  +  -r-  — 5- 

-,2  r  or     2  -,^2 

6v  r   dG 


Since  onlj'  complete  annular  or  circular  regions  are  to  be 
considered,  the  use  of  continuous  cyclical  functions  of  G  eliminates 
the  need  to  specify  boundary  conditions  along'  a  radial  boundary. 

In  the  present  work,  which  investip;atcs  the  buckling  of  annular 
sandwich  panels,  tlie  transverse  loading  function,  q,  is  not  considei-ed. 
Such  a  restriction  does  not  limit  the  application  of  the  obtained 
i-esiil'ts,  since,  for  snail  deflections,  the  ti-ansverse  loading  fiuiction 
does  not  influence  the  buckling  load  [3], 

Finally;  as  a  result  of  uniform  compression  along  cho  inner 
and  (oi-)  oi'.ter  boiindaries,  an  axisynunetric  buckling  mode  is  assvimed 
(B3)  to  result  from  the  lowest  critical  axial  load.   Such  an  assump- 
tion has  been  sIiov.ti  by  Olsson  [12]  to  be  valid  for  single-layer 
annular  panels. 

The  governing  equations  thus  become: 
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_  U    G   {A  [I   -f  (ru)]  f  t  V\}  -  ^  A  (,  N   ^)  =  0         (52) 

c  c  Lt   I  1'  <^'i'  J        J  r  dr    — rr  dr'^ 
c 

r  ?'      i,    cl'-'H  ^^i'^I  Ttl^u  1  du    u~l   „               ,^„^ 

c  L ^        cii'  I  ,   -^  L  ,  2  r  dr     2  j 

"-  c         -•  1-v     dr  r  -■' 


r.2   1  cl  ,   cK 

wnere      V   = r-  (r  -— ) 

r  dr    di- 


and , 

e 


-^  21;t  G-  j  -—+---+  X   —-^  6w  n   dQ  =  0  (54) 

<-'  L     c  c  Lt     2  di'_J   — rr  drJ      r 


^  illii  +  V   ^-  6u  n   dG  .-  0  (55) 

Ldr    1  rJ  r 


In  equations  (52-55)  the  subscript  r  has  been  omitted  from 
the  displacement  u,  witliout  confusion.   Also,  N   has  been  eliminated 
by  considering  Uie  pre-buckling  equilibrium  of  the  face  layers: 


dN  N    -  N.^ 
— rr  —rr   --9t    ^  /•cr-\ 
+ =  0                          (56) 

dr        r 


3.      A^riEil^  Stress   Distribution 

For    two-dimensional    axisymmotric   stress   distj'ibution    in  polar 
coordinates,    th.e  fjoveriiing    equations    are    [17]: 


d!"  .,   i.  ^  _   A  ^    0  (57) 

,    2        r    dr  2 

dr  r 


l^        ^   +    „      _  (58) 


-^•^'        1 


Vci7  +    ^f 
'f 
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2t  E    - 

f  f  /u      du\ 

^99  =  --7-  t  '    ''    d-rV  ^-^^) 


where  u  is  the  pre-bi!ckling  lateral  displacement. 
The  general  solution  of  equation  (57)  is: 


u  .=  C_r  +  C    /r  (60) 

i)     6 

where  C   and  C   are  arbitrary  constants  of  integration.   Substituting 

O  D 


(60)  into  (5S)  yields: 

N^.^.  -  --^  I  C  (l.v^)  -C   -^  (61.) 

l-v_^  -  r   -^ 

For  the  case  of  an  annular  panel  subjected  to  uniform  com- 
pression along  tlie  inner  and  outer  boundaries  (see  Figure  2),  boundai-y 
conditions  are: 


N    (a)  =  -  N. 
—yy  1 

(62) 

N      (b)   =:.   -  N 


where  a  and  b  arc  inner  and  outer  radii,  respectively,  and  N.  and 
N   are  inner  nr.d  outer  compressive  forces  per  unit  length,  respectively. 
Imposing  boundary  conditions  (62)  on  equation  (61),  we  obtain 
the  follO'VJng  axial  stress  distribut  j  01: 

N^r  =  "2  -  ^  <62) 


N^Q  .  -  -^  -  F  (64) 
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<v_     > 


Nc/^ 


Figure  2,   Annular  Sandwich  Panel 
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wliere  , 


and. 


D  . ^V-  ;   E  -.  -^ -i  (65) 


e  -  -  (66) 

b 


Since  K  _  and  X   are  considered  to  be  much  larger  than 

the  forces  produced  bj'  bonding  during  buckling,  the  axial  stress 

distribution  remains  essentially  unchanged  during  buckling. 

Also,  if   a  =  S  r;  0 ,  or   N  =  N  ,  it  is  easily  verified  that 
'  '       i     o  •^. 

equations  (63)  and  (64)  reduce  to: 


N   =  N..  =  -  N  (67) 

— rr   — 9G      o 


4.   Reduced  Equili'orium  Equations 

From  CHiaptcr  IV,  Section  2,  the  equilibrium  equations 
(reproduced  liere  for  convenience)  are: 

_  U   G   fl  [1  f  (ru)]  .  i   vU  -  i  f  (r  N    ^   .  0  (52) 
c  c  Lt  I  r  dr  J    r  dr  \   — rr  dr/ 

c 

c  Lt        dr    -.^1,2   r  dr    2 J 
c         -'   1-v.  "-dr  r  -■ 


(53) 


Zaid    [9]    and  Huang   and   Ebcioglu    [10]    unco\iple   siinilar    equa- 
tions  by   ope:'ating   on    tlie    second   witla   tlie   linear   opei-ator 

L(9)    =    -  -~    (r,o) 
T  r   dr        ' 
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and   eliminating    the   fii'st   large   coiiunon   bracket    in   both    equations. 
The  resulting   equation   is    then  directly    integrated   to    j-leld: 


1-v^  r'  r'  r' 


t.t    tE 
f    c      1 


—    Ir      N        -r—  drdr  +   -Tri'lnr-—     +— -r+   — 
r  ^    J    -rr   dr  2     \  2/2  r_ 

(68) 


where   C,  ,  C^ ,  and  C„  are  arbitrary  constants  of  integration. 
1 '   2 '       3 

Equation  (G8)  agrees  with  Huang  and  Ebcioglu' s  results  if 

Is    is  defined  by  equation  (67)  instead  of  (63).   It  sliould  be  noted, 
— rr 

hov.-ever,  tliat  tJie  uncoupling  procedure  described  above  yields  five 
constants  of  integration,  whjle  only  four  boundarj'  conditions  are 
available.   This  is  a  dii'ect  result  of  the  uncoupling  procedure  wliicli 
initially  increases  tiic  order  of  equation  (53). 

As  a  consequence  of  the   above-mentioned  complications  and 
other  considerations  '.vhich  will  be  discussed  later,  a  different 
technique  is  employed  in  the  present  work  to  uncouple  equations  (52) 
and  (53), 

Multiplying  equation  (52)  by  r  and  integrating  directly  without 
the  aid  of  eq\iation  (53)  ,  we  obtain: 


r~t   t    N   -,      C^   . 

c   — rr  I       11  .,,^. 

I-  ^     2tG  -^     2?G   ' 


where  ^ 

dw  ,    .^ 

dr 

and  C   is  an  a^bi  tvary  constant  of  iiitegration. 


37 


Comparing  equation  (69)  with  equation  (68),  the  advantage  of 
simplicity  becomes  inmiediately  evident.   Also,  it  will  be  made  clear 
in  the  following  sections  that  the  present  procedure  facilitates  the 
application  of  boundary  conditions,  and  suggests  analogies  between 
the  classical  single-laj^er  plate  theory  and  the  present  analysis  of 
sandwich  panels. 

Substituting  equation  (69)  into  equation  (53),  with  the  aid  of 
(63)  and  (70) ,  we  obtain,  after  some  simplification: 

3 

2  BC  T 

(,^0)r2  ^  ,    (,2_3c)  r  ^  ,    (jlr^  -  Ir^  .  3G)cp  .  -  -i-        (71) 
dr 


where,. 


c.JI,  -f  .  -^S!.  r.A-l 


t  £^G         G  (l-v^) 
c    c  c     f 

f  c  f 


and  D  and  E  are  defined  by  equation  (65). 

Equations  (60)  and  (71)  represent  the  uncoupled  equilibrium 
equations.   Since  (71)  is  a  second  order  differential  equation,  two 
constants  of  integration  are  generated.   Together  with  C   and  the 
constant  introduced  through  the  integration  of  (70) ,  we  have  four 
arbitrary  constants  of  integration  to  be  determined  bj^  boundary  con- 
ditions at  the  iuiier  and  outer  edges  of  the  annular  panel. 
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5.       Bounclai-y   Conditions 

From   Cliapter   IV,    Section   2,    the  boundary    conditions    (reproduced 
here   for   convenience)    are: 

/{2U   G     [f  +   I  ^"j  +   N        ^}   5w  n      de   .   0      (54) 
^     l        c    c  Lt  2    dr J        — rr    drJ  r 

c  -' 

!  {B  +  ^f  ?}  '^  \  ^'^  =  0  <^5) 

Integral  (51)  requires  the  specification  of  either  the  trans- 
verse deflection  w,  or  tlie  resultant  shear  stress  along  the  inner  and 
outer  edges.   Similarly,  integral  (55)  requires  the  specification  of 
either  llie  relative  lateral  movement  of  one  face  witli  respect  to  the 
other,  or  the  moment  produced  by  tensile  stresses  on  one  face  and 
compressive  stresses  on  t))o  other,  along  the  inner  and  outer  edges. 
Since  we  are  neglecting  the  bending  rigidity  of  the  faces,  the  large 
bracket  in  integral  (55)  represents  the  total  edge  moment,  wliile  u 
becomes  analogous  to  the  slope,  cp,  used  in  the  formulation  of  boundary 
conditions  for  tlie  classical  theory  of  single-laj^er  circular  plates  [3], 

Specifically,  in  the  present  work,  integrals  (54)  and  (55) 
are  satisfied  thi'oug]i  t'.io   following  choice  of  boundarj'  conditions  at 
the  inner  and  outei'  edges: 


ptt    N  ,  -, 
At   r  =  a:   u  +  j  -— -  +  -^^  qi  =  0  ;     u  =  0  (74) 

•-  ^     2tG  -^ 


At   r  =  b:        w  ;.:  0  ;     u  =  0  (75) 
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Conditions  (75)  are  analogous  to  boundary  conditions  termed 
clamped  or  built-in  in  the  classical  thcoi'y,  while  conditions  (74), 
termed  "slidei'"  in  tlie  present  v/ork,  have  previouslj'  been  emploj^ed 
for  stabilitj-  analysis  of  sandwich  columns  [1]  and  single— layer 
annular  panels  [12].   Physically,  such  a  restriction  could  be  approx- 
imated by  allov/ing  a  shaft  or  rigid  cylinder  to  occupy  the  central 
hole  (see  Figure  2). 

As  tlie  inner  radius,  a,  shrinks  to  zero,  conditions  (74) 
become  identical  to  boundary  conditions  present  at  the  center  of 
a  circular  san.dv.ich  panel,  witliout  central  hole,  constrained  alo)"»g 
the  outer  edge  only.   This  limiting  process  provides  a  check  for  our 
final  results  since  the  stability  problem  associated  with  such  a 
panel  yields  a  relatively  simple  solution. 

6.   Comparison  witli  Other  Theories 

Coir.parison  of  the  present  work  with  existing  theories  can  be 
facilitated  through  the  use  of  Table  1.   Referring  to  this  table,  an 
analogy  between  single-layer  theory  and  sandwich  theory  is  observed. 
Meissnei-' s  [3  1]  equation,  wliicl:  can  be  solved  by  means  of  bessel 
functions,  is  leduced  to  a  homogeneous  differential  equation  (1:  =  0) 
if  the  sliear  I'esultant  is  made  to  vanish  at  any  arbitrary  radius. 
This  becomes  immediately  evident  if  we  compare  the  boundary  condition 
which  specifies  zero  shear  with  the  second  form  of  the  equilibrium 
equation.   Consequent Ij",  since  the  shear  at  the  center  of  a  circular 
panel  must  varnish  because  of  symmetry,  Timoshenko'  s  [3]  bessel  equa- 
tion of  oi'dei  one  is  a  horaogeneous  differential  ccuation. 
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Analogously,  one  of  two  equilibrium  equations  describing 
an  annular  region  for  a  sandwich  panel  reduces  to  a  homogeneous 
equation  (C  =  0)  if  tlie  shear  resultant  is  specified  to  vanisli  at  some 
arbitrary  radius.   Such  a  simplification  also  provides  a  direct 
correspondence  between  u  and  cp  (see  equation  (69)).   In  the  case  of 
a  circular  sandwicli  panel,  a  homogeneous  equation  again  results  as 
a  consequence  of  the  symmetry  involved,  and  one  of  the  two  equilibrium 
equations  yields  bessel  functions  of  order  one,  as  is  the  case  in  the 
single  layer  theoiy. 

The  second  equilibrium  equation  for  a  circvilar  sandwich  panel, 
which  is  attributed  to  Huang  and  Ebciog]n,  is  deduced  from  equation  (69) 
ratlier  than  the  original  form  in  which  it  appeared  (see  equation  (68)). 
Consequently,  without  employing  the  present  uncoupling  technique,  the 
above  stated  analogies  would  not  be  evident. 

Since  the  present  work  deals  with  an  annular  sandwich  panel, 

it  becomes  obvious  tliat  bessel  functions  cannot  be  employed  unless 

N  =N   (see  Equations  (6r;)  and  (67)).   Solutions  for  this  special  case 
o    i 

and  the  more  general  case  are  olitained  in  Chapter  V. 


aiAPTER  V 
AXI SYMMETRIC  BUCKLING  OF  AXNTJLAli  SAJvTiWICH  PANELS 

1.   Uniform  Axial  Sti'ess  Distribution 

If  the  pressure  along  tlie  inner  edge  of  the  amuilar  sandwich 
panel  is  equal  to  the  pressure  along  the  outer  edge,  or: 

N.  ==  N  (76) 

1     o 

then,  from  (65),  (67)  and  (72): 


D:rG=0;   1  =  1;   E=N 

o 

N  BN 

^rr  —  ^^o  '   ^^  -  ^  -  -F  '   "  =  Ar(NT2y 

o 

wh.ich  reduces  equations  (71)  and  (69)  to: 


(77) 


o  ^2.         ,        r-     BN  on                  B<^i^ 

2   d  cp        av                     o  2            I                           1 

^    2         dr        La-(N   /2)  J^              A-(N   /2) 

dr                                  o  o 

pft             N     -,  C 

c             o  11 

"  =^    "      "^ 7~  ^'   '    ~~  r 


(78) 


(79) 


Applying   the   first    of   boundary    conditions    (7-1)    to   equa- 
tion   (79)    immediately    reduces    (78)    and    (79)    to: 


2  r-   HN    r^ 

2    d   cr,  dp        To  ,    ,  .,         .  (80) 


■  •^  ,1-     l-iiN     1-  -I 

Ti  -^  ^-  d7  ^  La::Tn-;^-  "  'J  ^'  ^  ' 

dr  o 
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^■H  N  -j 

u  .  -I-/--  — °-  cp  (81) 

■-       2tG  - 
c 

Since  (80)  is  Bessel's  equation  of  ordei'  one,  a  general 
solution  of  equations  (SO)  and  (81)  takes  the  follov/ing  form  [IS]: 

1  1_ 

r-         RN  ,2      -.  r-      BN  2      -, 

^-wliirT^)  d ^ Vi [(a^tIt^)  ^J  '''^ 

1  1 

p{;t  ^     -]    r  r/      BN  .2  -,  r/      BN  x2  -I  -, 

•■  -  -  [-^ -  — ]  k\  [(sttttIt)  ■■] * \h  [(^woTTIt)  ']  }        «-^' 


2tG 


where  J   and  Y   are  bessel  functions  of  oixler  one  of  the  first  and 
second  kind,  respectively,  and  A   and  A   are  arbitrary  constants 
of  integrat  i  on , 

If  we  )iov."  impose  the  remainder  of  boinadai-y  conditions  (74) 
and  (75)  on  equations  (82)  and  (83),  we  have: 


-tt  <K    ) 


-[-2-"lt-][Vl<^^'>^-    Vl<^S)]=0 


r-U  (N    )        -)  |--  -, 

-h~'---T;;-"jL-Vi^^^- Vi(^'M=o 


(84) 


2tG 
whei'e 


.B(K   )  ,2 

'         o   cr 


^  =    [A~iOfT"/2l)      ^    ;       ^    -    b  ^«^^) 


o   cr 
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These  equations  can  be  satisfied  by  taking  A  =  A  =0.   Then  the 
deflection  at  each  point  of  the  panel  is  zero  and  we  obtain  the 
trivial,  undeflectcd  form  of  equilibrium  of  the  panel.   The  buckling 
form  of  equilibrium  of  the  panel  becomes  possible  only  if  equations 
(84)  yield  values  for  A   and  A   different  from  zero,  which  requires 
that  the  determinant  of  tlie  coefficients  of  these  constants  vanish. 
Therefore,  after  multiplying  bj'  (4/t  t)  and  considering  the  first 
of  equations  (73),  our  critical  condition  becomes: 


[(x)   -  2]   [j^(,OY^(^P)  -  J^  (^P)Y^(^)J=  0 


(86) 


Equal  ion  (86)  closely  resembles  tlie  critical  condition  obtained  by 
Olsson  [12]  for  a  single-layer  panel  constrained  in  a  similar  mannei' 

For  a  given  value  of  3,  the  smallest  corresponding  value  of 
[1  for  wliicli  the  second  large  bracket  in  equation  (86)  vanishes  is 
given  in  Table  2  below  [12]. 


TABLE  2 
LOV/EST  V.ALUE  OF  [i.   SATISFYING  EQUATION  (86) 


li 

0.0000 
3.  832 

0.0256 
3.840 

0.0526 
3.860 

0.0909 
3.924 

0.1000 
3.  942 

p 

M- 

0.1111 
3.96G 

0.1250 
4.000 

0.1433 
4.045 

0.1667 
4.116 

0.2000 
4.235 

0.2500 
4,44  5 

0.3333 
4.905 

0.3956 
5.355 

0.5000 
6.39  1 

0.5461 
7.016 

l-i 

0. G2S5 
8.523 

0.6897 
10.175 

0.7634 
13.312 

0.8333 
18.87  3 

1.000 

CO 
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And,    from  the  first   of    (85): 


N  .  2 

o\  u. 


(f 


A/cr        ^2^      2.„ 
b   B+p,  /2 


(87) 


It  can  be  shown  that  (N  /A)    given  b\'  equation  (87)  is  a  monotonic 
o   cr 

increasing  function  of  ji.   Therefore,  the  lowest  value  of  p.  results 

2 

in  the  lowest  value  of  (N  /'A)   .   Also,  since  b  B  is  always  positive, 

o    cr 

equation  (87)  yields  values  of  (X  /A)    which  are  greater  than  or 

o    cr 

equal  to  two;  and  we  can  therefore  conclude  that  tlie  lowest  root  of 
equation  (86)  is  always  given  by  equation  (87)  in  conjunction  with 
Table  2.   (See  Figure  3.) 

The  first  of  boundary  conditions  (75)  has  not  been  used 
to  obtain  tlic  above  I'esults  since  the  buckling  load  is  independent  of 
a  transvei-se  translation  of  the  panel  as  a  whole.   In  this  respect  our 
analysis  parallels  the  classical  single-layer  theory  [3], 

Foi'   a  ^  p  =  0,  equation  (87)  becomes: 

cr    b  B+(3.832)  /2 

Equation  (88)  agrees  with  Huang  and  Ebcioglu' s  [lO]  results  if  the 
present  notation  is  used.   The  validity  of  this  lim.iting  process  was 
discussed  in  Section  5  of  Chapter  IV. 

As  G   approaches  a  vei'y  large  value,  tlie  first  term  in  the 
denominator  of  equation  (87)  becomes  dominant,  and  we  obtain: 


b 
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where  the  "effective  bending  rigidity,"  d' ,  is  defined: 


£2  2   „ 

t  t  t^E^ 
T.'  off 
D  =  5—  (90) 

2(l-v_^) 


Equation  (89)  in  conjimction  with  Table  2  agrees  with  the 
results  obtai)icd  by  Olsson  [12]  for  a  single-layer  annular  panel 
subjected  to  uniform  inner  and  outer  axial  pressures  of  equal  intensity 

For   a  =  p  =  0,  equation  (89)  becomes: 


,„  X      (3,832)^d' 

(\^r  =  ■ 2 (91) 

b 


Equation  (91)  is  identical  to  the  buckling  load  obtained  by 
Timoshonlco  [3]  for  a  circular  single-layer  panel  clamped  along  the 
outer  edge. 

This  limiting  process  is  intuitively  exi^ectcd,  since, 
for  G  =  0°,  the  only  structural  function  of  the  core  is  to  control 
the  distance  between  the  face  layers.   A  similar  relationship  exists 
between  the  web  and  f  ]  aiigcs  of  an  I-beam. 

2.   Genera]  Solution 

We  nov/  return  to  a  general  solution  of  equations  (69)  and 
(71).   Inl  reducing  the  djmensionless  variable,  T],  through  the  trans- 
formation: 

Tl  =   r/b  (92) 
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equation  (71)  becomes: 

2         ^  d\  2     d-o       4    2  ^^^^1^^ 

T1^N11"+1)  ?-|  +  TlCNll'^-S)  2^  +  (RT]^_Qn^3)cp  =  -  ~^-         (93) 


where, 


2 
N  =  --   ;       R  =  b  HX  ;      Q  =  IN  (94) 


SiKilaxly,  equation  (69)  becomes: 


-^t     N   -I       C 

C     —IT   I  1 


ptt     N    -I        C    ^ 


Since  ue  seek  a  scries  solution  of  equation  (93),  and  choose 
to  expand  our  series  about  the  point  T|  =;  1 ,  the  following  additional 
transf ormatioii  is  intrcxluced: 

§  =  'q  -  1  (9G) 

The  reasons  for  seeking  a  solution  about  the  point  T|  -  1  (§  =  0) 
will  be  discussed  in  Sections  5  and  6  of  tliis  chaptei'. 

V/i  th  tlie  aid  of  equation  (96),  and  after  some  simplification, 
equation  (93)  becomes: 

[(N^l)  +  2(2X+1)^  +  (CN+l)?;^  +  4N§^  +  Ng^]  ^  +  [  (N-3) 

+  3(K-1)§  +  3N§^  +  N§^]$  +  [(R-Qf3)  +  (4R  -  2Q)§ 

2       3      4        bV  (1+5) 
+  (6R--Q)§   +  4R5   +  R§  ]cp  zr ^ (97) 
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To  obtain  a  complement ai-y  solution  of  equation  (97)  the 
following  infinite  series  is  employed  [IS]: 

CO 

9^  =   ^  A  E,^  (98) 

"^    k=0   '^ 

where  the  coefficients  A   arc  functions  of  the  elastic  and  geometric 
properties  of  the  panel   and  the  critical  buckling  load.   The  radius 
of  convergence  of  series  (98)  can  be  shown  (see  Section  6,  Qiapter  v) 
to  be  of  sufficient  magnitude  for  our  particular  problem. 

Substituting  (98)  into  a  homogeneous  form  of  equation  (97), 
and  collecting  coefficients  of  common  powers  of  ^j  ^'e  have: 

{[2(N+1)A2  +  (N-3)A^  +  (R-qH3)A^]  +  [G(Nil)A^  +  2i5:;-l)A^ 

+    (R-Q+3X)A,  +  (4R-2Q)A  ]§  +  [12(N+1)A^  +  (27Xf3)A^ 
1  o  4  3 

4  (R-Q+}8X-1)A„  +  (4R-2Q+3N)A   +  (6R-Q)A  ]t^ 
2  1  o 

+  [20(X+1)A^  4  4(13N+3)A^  +  (45XfR-Q)A^ 
5  4  3 

+    (14X+4R-2Q)A„  +  (X+6R-Q)A   +  4RA  ]§ 
2  1       o 


}   - 


(99) 


In  order  that  this  series  vanish  for  all  values  of  §  in  some  region 
surrounding  §  =  0,  it  is  necessary  and  sufficient  that  the  coefficients 
of  each  power  of  §  vanish  [18].   This  pi'oduces  the  following  relations 
in  which  some  coefficients  have  been  eliminated  tlxrough  tJie  accumulative 
introduction  of  pieviovisly  computed  coefficients: 
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(R-Q+3)     .  (N-3) 

r(5N-l)  (R-Q+3)  (4R-2Q)'l  r(5N-l)(N-3)       (R-Q+SN)"!   ,  ,^  _^  , 

A      -  f-   IHzSl-         (R-Q+3)  (R-Q+lSy-l)  (27N+3)(4R-2Q) 

^^L-12(N+1)    ^  24(N+1)^  '  72(N+1)^ 

(27N-;  3)  (5N-1)  (R-Q+3)~j  ^  ["_    (4R-2Q+3N)  (N-3)  (R-Q+lSN-1) 

72(N+1)^  J      o^L"         12(N+1)      '  24(N+1)2 

^    (27N^3)(R-Q+3X)    _    (27N+3)  (5N-1)  (N-3)"1  ^  ^^^^^ 

72(X+1)^  72(N+1)^  -*      """ 


r(13N+3)(6R- 
•-        60(N+1)^ 


-  I  "(13N+3)(6R-Q)  (13N+3) (R-Q+3) (R-Q+18N-1) 

^^L-         2         -  120(N+1)^ 


_    (13X+3)  (27.\+3)  (4R-2Q)  (13X+3)  (27X+3)  (5N-1 )  (R-Q+3) 

3  "^  4 

360 (X+1)  360 (N+1) 

(45X+R-Q) (5X-1) (R-Q+3)  (45X+R-Q) (4R-2Q)  (R-Q+3) (14X+4R-2Q) 

120(N+1)"^  120(X+1)^  40  (N+1) 

4R        ~[    ^  r(13X+3)(4R-2Qt-3X)    _    (13X+3)  (N-3)  (R-Q+18N-1) 

-20(N+1)J^^L  60(,.,,)2  -  i20(N+l)' 


(13N+3) (27N+3) (R-Q4  3N)  (13y+3) (27N+3) (5N-1) (N-3) 

3  ^  4 

360 (N+1)  360 (N+1) 

(45X+R-Q) (5N-1) (N-3)  (45N+R-Q) (R-Q+3N) 

120(N+1)^  120(N+1)^ 

14N+4R-2Q)  (N+6R-Q)"]    .  ,^_„, 

~~;2  20  (N+1)  J  \  ^^"•^'' 


(N-3) ( 

+   — 


40 (N+1)' 


\  = 


v.'here  A      and  A,    remain   arbitrary. 
o  1 
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Because  of  the  complexity  of  the  computations,  ancl  the 
immediate  requirements,  no  recursion  formula  is  sought  in  the 
present  analysis. 

With  the  aid  of  (100-103),  our  complementary  solution  becomes: 

rn  (F^         A  /i    (^--Q^^^)    f2  ^  r(5N-l)(R-Q4-3)    (4R-2Q)~]  ^3 

f-  (6^^-Q>     (R-Q+3)(R-Qfl8N-l)    (27N+3)  (4R-2Q) 
•'L"  12(N+1)  "^       o......^2       "^ 


2 
24(N+1)  72(iN+l) 


(27N+3)  (5N-1)  (R-Q+3)~l  ^4  ^  r(13N+3)(6 
72(N+1)        -•      ^   60(N+1) 


(6R-Q) 
2 


(13N+3) (R-Q+3) (R-Q+18N-1)  _  (13X+3) (27N+3) (4R-2Q) 
120(N+1)^  360  (N+l)*^ 

(13Ni-3)(27N-t3)  (5N-1)  (R-Q+3)    (45N+R-Q)  (5N-1)  (R-Q+3) 


4  3 

360 (N+1)  120(N+1) 

(45y+R-Q)  (4R-2Q)   (R-Q4  3)  (14X+4R-2Q)    4R ~|   5 

"    120(N+1)^     '       40(N+1)2       -20(N+1)_J^^ 

r     ^X-3)   _2   r(5N-l)(N-3)    (R-Q+3N)~j  p3 

r   (4R-2g:-3N)    (y-3)  (R-Q^,18N-1)    (27N+3)  (R-Q+3y) 
-"    12(N+1)   ^      ^^^^.^^^2      -     ^^^^^^^2 

(27N+3)(5N-1)  (N-3)~l  _4   r(13Ni-3)  (4R-2Q+3N) 
I  I      +  I 

72(N+1)       -      "-     60(X+1) 
(13N+3) (N-3) (R-Q+18N-1)    (13N+3) (27N+3) (R-Q+3N) 


3  3 

120(N+1)  360 (N+1) 

(13N+3)  (27.\+3)  (5N-1)  (N-3)    (45N+R-Q)  (5N-1)  (N-3) 

+ 4 3 

360  (N+1)  120(yVl) 

(45N+R-Q) (R-Q+3N)    (N-3) (14N+4R-2Q)    (N+6R-Q)"] 
^  120(^.1)2      '      40(N+1)2      "20(N;i)J 


(104) 
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Following  a  procedure  similar  to  that  eiaployed  in  the  previous 
section,  equations  (95)  and  (104),  together  with,  a  particular  solution 
of  equation  (97) ,  are  constrained  according  to  boundary  conditions 
(74)  and  (75). 

Imposing  tlie  first  of  boundary  conditions  (74)  on  equation 
(95)  eliminates  tlie  need  to  find  a  particular  solution  of  equation 
(97),  since  ve  have   C  =  0.   The  subscript  c  can  therefore  be 
eliminated  from  equation  (98)  and  our  solution  becomes  a  complete 
solution. 

From  equations  (92)  and  (96),  we  note  that  at  i'  =  b,  T]  =  1 

and  £■  =  0.   Therefore,  from  equation  (95)  ,  recalling  that  C  =  0 

and  N   (b)  =  -  N   (see  equation  (62)),  the  second  of  boundaj-y  condi- 
rr         o 

tions  (75)  becomes: 


[(f)  -  ^] 


cp  (§=0)  ■--  0  (105) 


And  since  (N  /A)  ?;^  2  in  general,  condition  (105)  imposed  on  equa- 
tion (104)  yields   A  =  0 
o 

Similarly,  the  second  of  boundary  conditions  (74),  evaluated 
at   r  i^  a  (^=3-1),  and  expressed  through  equation  (95)  becomes: 


[(— )  -  2lcp  (§  =  P-1)  =  0  (106) 


After  substituting  equation  (104)  into  (106)  with  A  =  0, 

o 

and  noting  that  in   general  (N./A)  ?^  2,  we  arrive  at  the  genei-al 
buckling  criteria  for  an  annular  sandwicli  panel: 
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.  /^^  ^^    (^^-3)   ..  ^.2   r(5N-l)(N-3)    (R-Q+3N)"|  .^    ^, 

r   (4R-2Q+3N)    (N-3) (R-Q+18N-1)    (27N+3) (R-Q+3N) 
'  L"   12(X-.l)    ^      24(K.1)2      '  "   72(N-M)^ 

(27N+3)  (5X-l)(N-3)~]^„  _4   r(13N+3)  (4R-2CH3N) 


3 
72(N+1) 


-^         •-     60(Nhl)' 


(13N+3) (N-3) (R-Q4I8X-I)    (13N+3) (27N+3) (R-Q+3N) 


120 (N+1)  360 (N+1)^ 

(13N+3) (27N+3) (5N-1) (N-3)  _  (45N+R-Q) (5N-1) (N-3) 
■*"                  4  3     " 

360 (N+1)  120(N+1) 


(45N+R-Q)  (R-CVh3N)    (N-3)  (14N+4R-2Q)    (N+6R-Q)  |..,  ^.5 
120  (N+1  )'^  40(N4l)^        ^^' 


N+6R-Q)~|  , 
0(N+1)  S'" 


1 
. I  =  0  (107) 


As  explained  in  the  previous  section,  the  imdeflected  form 
of  a  panel  in  equilibrium,  i.e.,  Wien  A  -  0,  is  of  little  interest. 
Therefore,  ajjproximate  critical  buckling  loads  of  the  panel  under 
aiscussion  may  be  computed  by  considering  a  finite  numbei'  of  terms 
in  equation  (107). 

The  first  of  boundary  conditions  (75),  as  stated  earlier,  is 
not  used  in  obtaining  the  above  results,  since  the  buckling  load  is 
independent  of  a  transverse  translation  of  the  panel  as  a  wliole. 


3.   Successive  Approximations 

From  equations  (94)  ,  (72)  and  (65) ,  we  have: 
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N  rr 


2A(l-32)     ^V^  ^i^ 


5^(N  -N.)    P^(N  -N.) 
o   1  o   1 


(108) 


R  = 


2^\      (V^^\) 


,2     (N  -N.) 
)        o   1 


(109) 


r(N  -X.)    3  (N  -N.) 
o   1  o   1 


(110) 


We  designate  the  first  approximation  to  be  that  whicli 
considers  only  the  first  term  of  the  infinite  series  (107),  and 
consecutively  add  a  term  foi'  each  succeeding  approximation. 
Therefore,  v.'itli  the  aid  of  (108-110)  ,  and  after  som.e  simplification, 
we  have: 

First  Approximation 

Tiie  first  approximation  yields  no  results  since  the  loadiiig 
functions  are  not  present  in  the  fii'st  term  of  series  (107). 

Second  Approximation 


4A(l-3") 


2[(N  )   -3  (N.)   ] 
o  cr      1  cr 


+  2 


3  [(N  )   -(K.)   ]     r[(N  )   -(N.)   ] 
o  cr    1  cr  o  cr    i  cr 


-  (p-1) 


2A(l-3^) 


[(N  )   -3"(N.)   ] 
o  cr     1  cr 


-  3 


le-'rdX  )   -(N.)   ]      3'^[(N  )   -(X.)   ] 

1-      o  cr    1  cr  o  cr   i  cr    -^ 


r-O  (111) 
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Third  Approximation 


2A(l-&  ) 


[(N  )   -fl^CN.)   ]    ""^ 


U^[(N  )   -(N.)   ]    3^[(N  )   -(N.)   ] 
i-     o  cr    1  cr         o  cr    i  cr 


2A(l-3  ) 


&  [(N  )   -(N.)   ] 
o  cr   1  cr 


[(N  )   -g^(N.)   ] 
o  cr     1  cr 

P^[(N  )   -(N.)   ] 
o  cr    1  cr 


+  1 


2A(l-3^) 


[(N  )   -P  (N.)   ] 
o  cr      1  cr 

P^[(N  )   -(N.)   ]    3^[(N  )   -(N.)   ] 
o  cr    1  cr         o  cr    i  cr 


(p-1) 


10A(l-p  ) 


5[(N  )   -r(N.)   ] 
o  cr     1  ci' 


P"[(X  )   -(N.)   ]    3"[(N  )   -(N.)   ] 
o  cr    1  cr         o  cr    i  cr 


2A(l-g  ) 


P  [(N  )   -(N.)   ] 
o  cr    1  cr 


[(N  )   -3^(X.)   ]      , 

-^L^£ i--^  -  3  I  (P-l)2 

3  [(N  )   -(N.)   ] 
o  cr    1  cr 


OA/n  r2^  [(N  )   -3  (N.)   ] 

2A(l-p  ) o  cr 1  cr 

P^[(N  )   -(N.)   ]  B^[(N  )   -(N.)   ] 
o  cr    1  cr        o  cr    i  cr 


+  1 


4A(l-g  J 


P^KX  )   -(N.)   ]    5^1  (N  )   -(N.)   ] 
o  cr    1  ci-         o  cr    i  cr 


2[(N  )   -P^(N.)   ]    2b^B[(N  )   -?^(N.)   ] 
o  cr     1  cr  o  cr     i  cr 

+ 


■[(N  )   -(X.)   ] 
o  cr    1  cr 


2         2 

-  2b  B   (5-1)   =  0 


(112) 


Fourth  Approximation 


72 


2A(1-S  ) 


[(N  )   -  P  (N.)   ] 
o  cr      1  cr 

?^t(N  )   -(N.)   ]     0^[(N  )^  -(N.)   1 
o  cr   1  cr  o  cr   1  cr 


-  3 


+  1 
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o A  f  T  _o  ^ 


^<^o>cr-g'<^i^r^    J 


P^[(N  )   -(X.)   ]    P^[(X  )   -(X.)   ] 
o  cr    1  cr         o  cr    i  cr 


2A(1-P  ) 


P  [(X  )   -(X.)   ] 
o  cr    1  cr 


[(N  )   -P^(X.)   ]     ""^ 
o  cr     1  cr 

H.  1 

P  [(N  )   -(X.)   ] 
o  cr    1  cr 


(8-1)  +  (12 


10A(l-3  ) 


LP  ^^\Kr-^\^c.^ 
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5[(N  )   -  B  (N.)   ] 
o  cr       1  cr 


(1-B^ 


[(N  )   -S  (N.)   ] 
o  cr     1  cr 


o  cr 1  cr 2A(1-R  )  o  cr i  cr 

[(K  )   -(N.)   ]  P^[(N  )   "(N.)   ]    B^[(N  )   -(N.)   ] 

■- ^  o  cr    1  cr-'       J  L'-'  ■■  \  Q-'cr    i  cr         o  cr   i  cr 


2A(l-3  ) 


[(N  )   -S  (N.)   ] 
o  cr     1  cr 

B^[(N  )^,.-(N.)^  ]    P^[(N  )   -(N.)   ] 
i-     o  ci'    1  cr         o  cr    i  cr 


+  1 


-  ]2 


2A(i-S  ) 


B  [(N  )   -(N.)   ] 
o  cr    1  cr 


[(N  )   -P  (N.)   ] 
o  cr      1  cr 

P^[(N  )   -(X.)   J 
o  cr    1  cr 


-I  2 


+  1 


4A(1-S  ) 


2[(N  )   -3  (N.)   ] 
o  cr     1  cr 


3  [(N  )  -  (N.)   ]    3  [(N  )   -(N.)   ] 
I-     o  cr    1  cr         o  cr    i  cr 


2b^B[(N  )   -3^(X.)   ] 

o  cr      1  cr     „,  2 

+ 2b  B 

P  l(X  )   -(X.)   ] 
o  cr         1  cr 


(S-1)"  + 


2A(l-p^ 


3  [(N  )   -(X.)   ] 
-      o  cr    1  cr 


[(X  )   -3  (X.)   ] 
o  cr     1  cr 

P^[(X  )   -(X.)   ] 
o  cr    1  cr 


-  2  r- 


+  1 


oAM  =2,  [(X  )    -r  (X.)    ] 

2A(l-3  )  o  cr      1  cr 


'[(X  )   -(X.)   ]    3  [(X  )   -(X.)   ] 
o  cr    1  cr  o  cr    i  cr 


e,  2^,  [(X  )   -3^(X.)   ] 

8b  B  o  cr     3  cr 

~Z2  Tex"!    ^~)    r 

p  o  cr    1  cr 


4b^B 


+  3 


2A(l-3.  ) 


o  cr    1  cr 


[(X  )   -p^(X.)   ] 
o  cr     1  cr 

3^[(X  )   -(X.)   ] 
o  cr    1  ci- 


+  1 


2A(l-3  ) 


[(X  )   -3  (X.)   ] 
o  cr     1  cr^ 


3  [(X  )   -(X.)   ]    3  [(X  )   -(X.)   ] 
-  -      o  cr    1  cr  o  cr    i  cr 


34A(l-[-.  )  o  cr      1  cr 


O'er  ^  -i'cr^    2b2B  ^  ^  V  cr"^  ^\>cr^ 


3^[(X  )   -(X.)   ]     3^[(X  )   -(N.)   ]      3^    f^^\^cr"^"'i^cr^ 
>-     o  cr    1  cr  o  cr    i  cr 


-  2b^B  -  1 


2A(l-3^) 


f(\^r-P  ^\Kr^ 


J   3^[(x  )  -(X.)  ]   3^[(x  )  -(X.)  ] 

-^   "-     o  cr    1  cr         o  cr    i  cj- 


+  1 


54A(l-3  ) 


27 [(X  )   -3  (X.)   ] 
o  cr  '^    1  cr 


+  3 


4A(l-3  ) 


■^  ^^Vcr-(^"i>cr3     P't(N^)^^,-(X.)^^J      -L^'^  (X^)^^,-(X.  )^^.] 
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2[(N  )   -S^CN.)   ]    2b^B[ (N  )   -3^(N.)   ] 
o  cr     1  cr  o  cr     i  cr 


[(N  )   -(N.)   ] 
o  cr    1  cr 


B  [(N  )   -(N.)   ] 
o  cr    1  cr 


2b  B   - 


54A(l-3  ) 


3^[(N  )    -  (N.)   ] 
o  cr     1  cr 


27  [(N  )  -?j    (N.)   ] 
o  cr      1  cr 


+  3 


10A(l-5  ) 


5[(N  )   -3^(N.)   ] 
o  cr     1  cr 


-  1 


[(N  )   -(N.)   ]         B  [(N  )   -(N.)   ]    P  [(N  )   -(N.)   ] 
o  cr    1  cr      J   -'^  o  cr    i  cr         o  cr   i  cr 


2A(1-B  ) 


[(N  )   -P  (N.)   ] 
o  cr     1  cr 

B^[(N  )   -(N.)   ]    B^[(N  )   -(X.)   ] 
o  cr    1  cr  o  cr    i  cr 


-  3  )(?-l)   =  0 


(113) 


Fifth  Approximation 


360 


2A(l-p") 


[(N  )   -B  (N.)   ] 
o  cr     1  cr 


+  1 


B  [(N  )   -(X.)   ]    p^[(N  )   -(N.)   ] 
o  cr    1  ci'         o  cr    i  ci' 


180 


r  2 

■  2A(1-B  ) 


[(N  )   -3  (N.)   ] 
o  cr     1  cr 


-  3 


4  1 


'[(N  )   -(N.)   ]    3"[(N  )   --(N.)   ] 
o  cr    1  cr         o  cr    i  cr 


2A(l-3^) 


8  [(i\  )   -(N.)   ] 
o  cr    1  cr 


[(N  )   -?  (N.)   ] 
o  cr      1  cr 

B^[(N  )   -(N.)   ] 
o  cr    1  cr 


(R-l)  -t  (60 


2A(l-p'') 


[(N  )   -(N.)   ] 
o  cr    1  cr 


[(N  )   -B'-CN.)   ] 
o  cr     1  cr 


-1  2 


+  1 


2  ■  -^  I 

-  [(N  )   -(N.)   ]     J 
o  cr    1  cr      -J 


10A(]-p  ) 


5[(N  )   -p'^CN.)   ]     1 
o  cr     1  cr      I 


P^[(N  )   -(In.)   ]    P^[(N  )   -(N.)   ] 
o  cr    1  cr         o  cr    i  cr 


2A(l-3  ) 


[(N  )   -p  (N.)   ] 
o  cr     1  cr 


P^[(N  )   -(N.)   ]    p^[(N  )   -(N.)   ]     J 
I-  o  cr    1  cr         o  cr    i  cr     -J 


60 


2A(l-p  ) 


.P  f<Vcr-^^\>cr^ 


[(N  )   -p^(N.)   ]     ]^  "■ 

o  cr     1  cr    ,1 
__ +  1  I 


P  [(X  )   -(X.)   ] 
o  ci-    1  cr 


4A(l-p  ) 


2[(N  )   -p  CN.)   ] 
o  cr    ^  1  cr 


L 


3"[(X  )   -(X.)   ]     P^[(N  )  „-(N.)_] 


o  cr    1  cr 


o  cr    1  cr 
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2b^B[(N  )  ,,-3^(N.)   ] 

o  cr     1  cr    „,  2 
+ 2b  B 

B  [(N  )   -(N.)   ] 
o  cr    1  cr 


(p-1)   + 


2A(l-p  ) 


'[(N  )   -(N.)   ] 
o  cr    1  cr 


[(N  )   -p  (N.)   ] 
o  cr     1  cr 

P^[(N  )   -(N.)   ] 
o  cr    1  cr 


+  1 


2A(l-3  ) 


[(N  )   -3  (N.)   ] 
o  cr 1  cr 

P  [(X  )   -(N.)   ]    p^[(N  )   -(N.)   ] 
>-     o  cr    1  cr     ^  •■  ^  o  cr    i  cr"^ 


8b^B[(N  )   -g^(lN.)   ] 

o  cr     1  cr     ^2 

+  2 ^^  ^ 

P  [(N  )   -(N.)   ] 
o  cr    1  cr 


+  15   - 


2A(l-3  ) 


P  [(N  )   -(N.)   ] 
o  cr    1  cr 


[(N  )   -  3^(N.)   ] 
o  cr 1  cr 

B^[(N  )   -(N.)   ] 
o  cr    1  cv 


+  1 


2A(l-?5  ) 


[(N  )   -(N.)   ] 
o  cr    1  cr 


[(N  )   -3  (N.)   ] 
o  cr 1  cr 

3^[(N  )   -(N.)   ] 
o  cr    1  cr 


-  3I  r  34A(i-3^) 


17[(N  )^^-3^(N.)   ]    2b^B[(N  )   -  3^(N.)   ] 
o  cr     1  cr  o  cr      1  cr 


'^^^'^c?c~.r^'^.">,J  3^[(N  )   -(N.)   ]       3^[(N  )   -(N.)   ] 

o  cr    1  cr  o  cr    1  cr  o  cr    i  cr 


2b  B 


+  5 


2A(l-3  ) 


[(N  )   -3  (N.)   ] 
o  cr      1  cr 


'  B"[(N  )   -(N.)   ]    p'^LCN  )   -(N.)   ] 
o  cr    1  cr         o  cy        1  cr 


54A(1-S  ) 


27 [(X  )   -3  (N.)   ] 
o  cr 1  cr 

^^f^^'o>r..-(^'>   ^     B^[(N  )   -(N.)   ] 
o  cr    1.  cr  o  cr    1  cr 


+  3 


4A(l-p  ) 


3  [(N  )   -(N.)   ] 
o  cr    1  cr 


2[(N  )   ■ 
o  cr 


MN.)   3    2b  B[(N  )   -B  (N.)   ] 

1  cr            o  cr    ^  i^cr-"      2 
4 2b  B 


3  [(N  )   -(.\.)   ] 
o  cr    1  cr 


3  [(N  )   -(N.)   ] 
o  cr    1  cr 


-  5 


2A(1-S  ) 


[(N  )   -3  (N.)   ]     ^ 
o  cr     1  cr    ,  ! 
—  +  1 


LP'[(^o),,-(N.)^^.]    b2[(N^)^^,-(N.)^^,]  ■  ':Lp'[(N,),,-(N,)^^] 


54A(l-3  ) 


27[(N  )   -?^(N.)   ]     If      ,,2 
o  cr __JL_l££._   3 10A(l-3  ) 


5[(N  )   -3  (N.)   ] 
o  cr      1  cr 


^  ^^^>cr-(^'i)cr^      -  -P  t^Vcr-^^^^r^     ^'"^^  ^  V  cr-^^^^cr^ 


-  1 
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2A(l-3  ) 


[(N  )   -3  (N.)   ] 
o  cr 1  cr 

3^[(N  )   -(N.)   ]    S^[(N  )   -(N.)   ] 
o  cr    1  cr         o  cr    i  or 


(P-1)" 


2A(1-S  ) 


[(N  )   -e  (N.)   ] 
o  cr      1  cr 


P^[(N  )   -(N.)   ]    3^[(N  )   -(N.)   ] 
>-      o  cr    1  cr  o  cr    i  cr 


+  1 


26A(l-3  ) 


[(N  )   -(N.)   ] 
o  cr    1  cr 


+  3  '  ' 

'[(N  )   -(N.)   ] 
o  cr    1  cr 


[(N  )   -3  (N.)   ] 
o  cr     1  cr 

^^^'^^r.K^-'^^-^.J        ^^L(N  )   -(N.)   ] 
o  cr    1  cr         o  cr   i  cr 


8b^B[(N  )   -p^(N.)   ] 

o  cr     1  cr     .,  2^  , 
+  5 4b  B  I  -  3 

3  [(N  )   -(N.)   J 
o  cr    1  ci- 


2A(l-3  ) 


[(N  )   -(N.)   ] 
o  cr    1  cr 


[(N  )   -3  (N.)   ] 
o  cr  1  cr 

3^[(X  )   -(X.)   ] 
o  cr    1  cr 


+  1 


26A(1-S  ) 


3  [(N  )   -(N.)   ] 
o  cr    1  cr 


13 [(N  )   -3  (N.)   ] 
o  cr     1  cr 

3^L(N  )   -(N.)   ] 
o  cr    1  cv 


+    3 


2A(]-3") 


[(N  )   -3  (N.)   ] 
o  ci'     1  cr 

B^[(N  )   -(N.)   ]    P^L(N  )   -(X.)   ] 
o  cr    1  cr         o  cr    i  cr 


34A(l-3^) 


3  [(X  )   -(X.)   ] 
o  cr    1  cr 


17 [(X  )   -p^(X.)   ]    2b^B[(X  )   -3^(X.)   ] 

o  cr      1  cr  o  cr      i  cr     „  2 

~:27t: — r—  ■"  ~2 2b  B  - 1 


[(X  )   "(X.)   ] 
o  ci-    1  cr 


3  [(X  )   -(X.)   ] 
o  cr    1  cr 


2A(1-3  ) 


[(X  )   -(X.)   ] 
o  cr    1  cr 


[(X  )  -3  (X.)  ] 

o  cr     1  cr 

3^[(X  )   -(X.)   ] 
o  cr    1  cr 


+  1 


26A(l-3  ) 


3  [(X  )   -(N.)   ] 
o  cr    1  cr 


13[(N  )  .-3  (X.)   ] 


L(X  )   -(X.)   ] 
o  cr    1  cr 


.al 


r  2 

54A(l-3  ) 

3^[(x  )  -(X.)  ] 

o  cr    1  cr 


27 [(X  )   -3  (X.)   ] 
o  ci"      1  cr 


'[(X  )   -(X.)   ] 
o  cr    1  cr 


+  3 


4A(l-32) 


3^L(X  )   -(X.)   ]    3^[(X  )   -(X.)   ] 
o  cr    1  cr  o   cr    i  cr 


2[(X  )   -B^CX.)   ]    2b2DL(X  )   -32(x.)   ] 
o  cr     1  cr  o  cr     i  cr 


3  L(X  )   -(X.)   ] 
o  cr    1  cr 
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2 

-  2b  B 


26A(l-p  ) 


13[(N  )   -3  (N.)   ] 
o  cr      1  cr 

P^[(N  )   -(N.)   ]     3^[(N  )   -(N.)   ] 
'-  o   cr        1  cr  o  cr    i  cr 


+  3 


r  2 

54A(l-p  ) 

P^[(N  )   -(N.)   ] 
o  cr    1  cr 


27 [(N  )   -3  (N.)   ] 
o  cr     1  cr 

P^[(N  )   -(N.)   ] 
o  cr    1  cr 


+  3 


10A(l-3  ) 


B  [(N  )   -(N.)   ] 
o  cr    1  cr 


5[(N  )   -?^(N.)   ] 
o  cr     1  cr 

P^[(N  )   -(N.)   ] 
o  cr    1  cr 


-    1 


2A(l-p  ) 


[(N  )   -3  (N.)   ] 
o  cr      3  cr 

;  I  p^[(N  )   -(N.)   ]    6^[(N  )   -(N.)   ] 
-I-     o  cr    1  cr         o  cr    i  cr 


=  ,.3 


2A(1-P  ) 


P  [(N  )   -(N.)   ] 
o  cr    1  cr 


[(N  )   -P  (N.)   ] 
o  cr     1  cr 

P^[(N  )   -(N.)   ] 
o  cr    1  cr 


■h  1 


88A(l-0  ) 


3^[(N  )   -(N.)   ] 
o  cr    1  cr 


44[(N  )   -S  (N.)   ]    2b  B[(N  )   -3  (N.)   ]      ^  ] 
o  cr      1  cr            o  cr      i  cr     „,  2  i 
— +  — 2b  B 


'[(N  )   -(N.)   ] 
o  cr    1  ci- 


'[(N  )   -(N.)   ] 
o  cr    1  cr 


10A(1--  ) 


3  [(N  )   -CN.)   ] 
o  cr    a  cr 


5[(N  )   -3^(N.)   ] 
o  cr     1  cr 

3^[(x  )  -(X.)  ] 

o  cr    1  cr 


2A(l-3  ) 


[(N  )   -3  "(X.)   ] 
o  cr     1  cr 


?^[(N  )   -(N.)   ]    3^[(N  )   -(N.)   ]     J 
1-     o  cr    1  cr         o  cr    i  cr     -" 


+  3 


2A(l-3.  ) 


[(N  )^-3  (N.)^  1 
o  cr    1  cr 


-  2 


3^[(N  )   -(N.)   ]    p^[(N  )   -(N.)   ] 
>-     o  cr    1  cr         o  cr    i  cr 


+  1 


88A(l-3  ) 


P"[(N  )   -(N.)   ] 
o  cr    1  cr 


44[(X  )   -3^(N.)   ]    2b^B[(N  )   -p^(N.)   ] 
o  cr     1  cr  o  cr     i  cr 


2b^B 


3  [(N  )   -(N.)   ] 
o  cr    1  cr 


3  [(N  )   -(X.)   ] 
o  cr    1  cr 


4A(l-3  ) 


3  '[(N  )   -(N.)   ] 
o  cr    1  cr 


2[(N  )   -3^(N.)   ]    2b^B[(N  )   -B^(N.)   ] 
o  cr     1  cr  o  cr     i  cr 

— -- — , +  _ 


'[(X  )   -(N.)   ] 
o  cr    X  cr 


2b^B 


[(X  )   -(X.)   ] 
o  cr    1  cr 


+  9 


2A(l-3^) 


[(N  )   -3^(X.)   ] 
o  cr     1  cr 


3  [(X  )  -(X.)  ]   3  [(X  )  -(X.)  ] 

o  cr    3  cr     ^  o   cr        i  cr 


2A(l-3  ) 


P.  [(X  )   -(N.)   ] 
•-     o  cr    1  cr 
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t(N  )   -3  (N.)^,J 
o  cr     1  cr 


-  3 


r 


24A(l-3  ) 


12[(N  )   -P  (N.)   ] 
o  cr     1  cr 


e^[(N  )   -(N.)   ]       I  P^t(N  )  ,,-(N.)_]     p^[(N  )   -(N.)   ] 
o  cr   1  cr     -  1-     o  cr   i  cr  o  cr   i  cr 


8b^B[(N  )   -S^(N.)   ]      _  ■ 
o  cr      1  cr     ,2. 

+  p 4b  B 

P  [(N  )   -(N.)   ] 
o  cr    1  cr 


18 


2A(l-3  ) 


3  [(N  )   -(N.)   ] 
o  cr    1  cr 


[(N  )   -3"(N.)   ] 
o  cr     1  cr 


-  -I-  1  ! 


P  [(N  )   -(N.)   ] 
o  cr    1  cr 


~^^12b^B[(N  )   -3^(N.)   ] 

^^£ L-££_  _  sb^B 

3  [(N  )   -(N.)   ] 
o  cr    1  cr 


^(3-1)^=0 


(11^) 

Computation  of  tlie  sixth  approximation  is  unnecessary  since 
it  can  be  shov.ii  that  the  fifth  apj^roximation  yields  acceptable 
results  (see  Section  5,  Cliapter  V). 

Equations  (111-111),  each  being  self-contained,  represent  the 
approximate  critical  buckling  criteria  for  an  annular  sandwich  panel 
constrained  by  bomidary  conditions  (74)  and  (75)   and  subjected  to 
uniform  radial  compressive  loads,  N.  and  N  ,  along  the  inner  and  outer 
edges,  i-especti  vely. 


4,   Numea'ical  Results  and  Discussions 

Once  tlie  inner  compressive  load,  N.  ,  is  pi-escribed  to  be  some 
multiple  of  the  outer  compressive  load,  N  ,  tlie  approximate  critical 

conditions  given  in  the  previous  section  are  comiiletely  defined  bj^ 

2 

three  dimcnsio)ilcss  parameters:   (N  /A)   ,  b  B,  and  p.   Thus  we  obtain 

o    cr 

algebraic  polynomials  in  (N  /A)    which  increase  in  degree  as  the  order 

o    cr 

of  apj)roximatiorx  increases  (ranging  from  first  degree  in  the  second 
approximation  to  fourth  degree  in  the  fifth  approximation).   For 
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obvious  reasons  we  consider  only  the  lowest  positive  value  of  (N  /A) 

o    cr 

satisfying  eacl^i  polynomial. 

Three  possible  loading  conditions  are  analyzed  in  the  present 

work:   (1)  N  =  N.  ,  (2)  N  =  0,  and  (3)  N  =  0.   However,  it  sliould  be 
o     1         o  1 

noted  that  the  techniques  employed  in  these  examples  are  applicable 

foi'  any  ratio  of  X   to  N.  . 
ox 

N  =  N. 


If  the  inner  and  outer  axial  compressive  loads  are  equal ,  then 
the  exact  solution  is  given  by  equation  (87)  in  conjunction  with 
Table  2,   Thus,  tlie  approximation  techniques  employed  in  the  previous 

section  ai's  unnecessai-y.   For  this  case.  Figure  3  shows  the  existing 

2 

relation  betweeji  (X  ,''A)    and  3  for  vai'ious  values  of  b  B, 
o    cr 


If  the  inner  edge  alone  is  subjected  to  axial  compression, 

then,  since  equations  (111-114)  yield  no  positive  values  of  (N  /A)   , 

o        cr 

it    can  be   concluded    that   buckling   never   occurs.      This   wo\ild   seem 
reasonable,    since,    from   eq^iations    (63-65),    such   a   reduction   I'esults 

in   a   relativelj'   large   tensile  N        compared  v.ith   a   relatively   small 

— y  B 

compressive  N_  .   Analogously,  for  a  rectangular  single-layer  panel 
subjected  to  compression  along  opposite  edges  and  tension  along  adja- 
cent edges,  Timoshenko  [3]  has  shown  that  a  large  tensile  load  will 
prevent  a  significantly  smaller  compressive  load  from  causing  insta- 
bility.  Ho'vevev,  '.\'e  must  keep  in  luind  that,  while  N   anrl  X   can 
^  '  '  '        -XX     -yy 

be  varied  j.ndependeni  ly  in  a  rectangular  panel,  such  is  not  the  case 
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for  a  circular  or  annular  panel,  since  the  following  equilibrium 
condition  must  be  maintained: 

dN     N   -  N„^ 

_i2Z  +  :2Z =90  ^  0   (56) 

dr        r 

N.  =  0 

1 

If  onlj'  the  outer  edge  is  subjected  to  axial  compression, 

the  second  approximation,  equation  (111),  can  be  solved  exiplicitly 

for  (N  /A)   : 
o    cr 


(t) 


^"^°'"-"  (115) 

cr    3+23  +  35 


2 

Since  b  B  does  not  enter  into  relation  (115),  it  is  obvioiis 

that  furtlier  approximations  must  be  considered. 

Due  to  the  complexity  of  the  calculations  involved,  a  graphical 
solution  is  employed  for  the  succeeding  approximations.   In  Figure  4 

approximations  two  througli  five  are  compared  bj'  plotting  (N  /A) 

2 
versus  3  foi-  various  values  of  b  B. 


All  quantities  appearing  in  Figures  3  and  4  are  dimensionlcss. 
Buckling  loads  may  be  obtained  in  the  appropriate  dimensions  by  using 
tlie  first  of  relations  (73). 
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Lo 


Figure  3.   Minimum  Critical  Values  of  (N  /A)  for  N  =  N 

o         o    i 
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^'approx.  Ci'iy 

r 


'^T^f+Vpprox;.   6^e>^  /^^^ 


T r 


p--Vh 


-T" 

.8 


,^         /.o 


Fijxure   4.      Minimum   Critical    Values   of    (N  /A)    for  N.    =    0 
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The  following  conclusions  may  be  drawn  from  Figures  3  and  4: 

1.  All    critical  values  of  (N  /A)  are  less  than  or  equal  to  two. 

o 

(Figures  3  and  4) 

2 

2.  As  the  value  of  b  B  approaches  infinity,  the  results  approach 

those  obtained  for  a  single-layer  panel  in  equation  (89). 
(Figure  3) 

3.  For  6  -   0,  the  present  theory  coincides  with  Huang  and  Ebcioglu's 
results  (equation  (88))  for  a  circular  sandwich  panel.   (Figures  3 
and  4) 

4.  The  second  approximation  in  Figure  4  is  the  exact  solution  for 
b^B  =  0. 

5.  An  annular  sandwich  panel  subjected  to  axial  compression  along  the 
outer  boundary  becomes  stronger  if  an  equal  compressive  load  is 
also  applied  along  the  inner  edge.   (Figures  3  and  4) 

6.  In  Figure  4,  the  third  approximations  yield  more  accurate  results 
than  the  fourth  approximations.   This  peculiarity   and  the  error 
bound  associated  with  the  fift]i  approximations  will  be  discussed 
in  tlie  next  section. 

7.  A  dual  response  is  apparent  in  Figure  4.   As  the  hole  increases 

in  I'olative  size,  the  panel  may  become  weaker  or  stronger  depend- 

2 

ing  on  tlie  value  of  b  B  and  the  range  of  P  being  considered. 

Such  a  beliavior  is  possible  because  both  tlie  shear  a)id  the  bending 
stiffness  of  a  sandwich  panel  enter  into  the  analysis.   An  annular 
single -layer  panel,  which  can  be  described  by  only  two  dimension- 
less  parameters,  exliibits  no  sucli  dual  response. 
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8.   If  buckling  loads  for  values  of  b  B  or  ratios  of  N   to  N  ,  not 

o      1 

considei'G'd  in  Figvires  3  and  4,  are  required,  equation  (]  14)  may 
be  used  directly.   However,  when  considering  only  values  of  p 
gi'eater  than  one-half,  the  third  approximation,  equation  (112), 
yields  acceptable  results.   (Figure  4) 

The  above  conclusions  are  valid  only  for  tlie  special  case  in 
wliich  boundary  conditions  (74)  and  (75)  are  applied. 

5.   Error  Refund 

Since  v/c  chose  to  expand  tlie  solution  of  equation  (93)  about 
the  point  T|  n.  1 ,  the  speed  of  convergence  of  the  series  solvition 
obtained,  equation  (104),  depends  on  the  proximity  of  the  entire 
annular  region  to  that  point.   Clearly,  as  the  hole  increases  in  size 
(p  approaches  one),  the  solution  converges  more  rapidly.   This  fact  is 
evident  from  Figure  4,  and  also  from  equation  (107).   Indeed,  it  can 
be  concluded  that  the  speed  of  convergence  is  the  slowest  v.lien   3^0. 

From  equations  (111-114)  the  approximate  critical  buckling 
parameter,  (N  /A)  _,  can  be  solved  exi^licitly  for  tlie  degenerate  casf; 
of   3=0.   The  results  of  this  simplification  are  found  in  Table  3. 

Howe^'er,  for  this  special  case,  the  exact  solution  (reproduced 
here  for  convenience)  is  available  from  Section  1  of  tliis  chapter: 

\\ 1^-684         ^ 

A  J         "      2 

cr    b  B  +  7.342 


TABLE    3 


APPnoXir.iATE   VALUES   OF    (N  /A) 

o    cr 


FOR  B  =  0 
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Order  cf 
Approximation 


(N  /A) 
o    cr 


First 

(No  results) 

Second,  equation  (111) 

2 

Third,  equation  (112) 

12 

b^B  +  6 

Fourth,  equation  (113) 

10 

b^B  +  5 

Filth,  equation  (114) 

15 

b^B  +7,5 

If  we  conpare  equation  (88)  with  Table  3,  it  can  be  concluded 

that  the  fifth  approximation,  for  3=0,  is  witliin  2.2  per  cent  of  the 

2 

exact  solution,  even  for  large  values  of  b  B.   And,  since  the  series 

solution  converges  more  rapidly  for  other  values  of  3,  equation  (114) 

yields  results  that  lie  within  2.2  per  cent  of  the  exact  solution  for 

2 

all  valties  of  b  B  and  p. 

Following  the  same  reasoning  outlined  above,  it  can  be  concluded 
that  the  third  approximation  yields  more  accurate  results  than  the 
fourth  approximation.   This  peculiarity  can  be  easily  verified  by 
ex-ijanding  a  series  solution  of  equation  (78)  about  the  point   r  =  b 
(5  -  0).   If  this  is  done,  the  resulting  approximations,  as  expected, 
coincide  exactly  with  those  listed  in  Table  3. 
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6 ,   Remarks  on  Convergence 

The  series  solution  of  an  ordinary  differential  equation 

possesses  a  radius  of  convergence  at  least  as  great  as  the  distance 

from  the  point  of  expansion  to  the  nearest  singularity  [19]. 

Equation  (93)  possesses  two  regular  singularities,  one  at 

2 

T|  =  0,  and  another  at  (nT]  +1)  =  0;  and  its  solution  was  expanded  about 

the  point   T]  =  1.   Therefore,  it  must  be  demonstrated  that  these 
singularities  do  not  inhibit  the  validity  of  our  solution  throughout 
the  entire  annul ai-  region  of  the  panel. 

Referring  to  Figure  5,  it  becomes  obvious  that  the  singularity 
at  T|  =  0   does  not  restrict  tlie  required  radius  of  convergence,  regard- 
less of  the  value  of  p.   It  is  therefore  necessary  only  to  show  that 

2 

the  singularity  occurring  at  (>n]  +1)  =  0,  lies  outside  the  annular 

region  and  its  reflection  illustrated  in  Figure  5,  for  all  values  of  p. 


-P 


^  ^e-fleciecl  I  maae    <}fcii\n.u let r 

^    region    -Prom    h-^1 


! T ^  H 

P  1  2-/3  2  ' 


Figure  5.   Radius  of  Convergence 


2 

With  the  aid  of  equation  (108),  (N1]  +1)  =  0  becomes: 


2 

^       [2A(1-S^)  -  (N  -P^N.)]  +1=0  (IIG) 


o       1 
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Therefore,  the  position  of  the  singularity  asfociated  with  equation 
(116)  depends  on  the  value  of  the  critical  buckling  load  which  becomes 
known  only  after  the  solution  is  obtained. 

For  tlie  particular  case  in  which   N.  =0,  (116)  is  satisfied 
if 

UJcr    ^2_p2 

However,  from  equation  (115)  and  Figure  4,  it  is  apparent 
that,  for  N.  =  0, 

-'  ^      2(l  +  R)(3-«) 


s  "^"""^  ^  ''  (118) 

A  /cr  2 

3  +  23  +  3p 

for  all  values  of  b  B, 

In  oi-der  tliat  conditions  (117)  and  (118)  be  satisfied  simul- 
taneously, Ti   must  satisfy  the  following  inequality: 

^2  ^  , 2(3^,^)3! (,,,^ 

[2(3-0)  -  2(1-3) (3+2P+33  )] 

2  2   • 

Equation  (119)  constrains  T]   to  be  greater  than  (2-3)   for 

all  values  of  3  ranging  from  0  to  1.   It  can  therefore  be  concluded 
that  the  singularity  associated  with  equation  (117)  lies  outside  the 
annular  region  and  its  reflection  illustrated  jn  Figure  5.   Similarly, 
it  can  also  be  shov.ii  that,  for  N  =0,  the  singularity  associated  with 
equation  (116)  lies  outside  this  ci^itical  region. 
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Thus,  the  radius  of  convergence,  p,  of  the  series  solution  of 
equation  (93),  expanded  about  the  point  T]  =  1,  is: 

p  ^  1  -  P  (120) 

which,  as  illustrated  in  Figure  5,  is  large  enough  to  encompass  the 

entire  annular  region  of  the  panel. 

Care  must  be  taken,  however,  when  imposing  boundary  conditions 

other  than  (74)  and  (75)  on  the  solution  of  equation  (93).   Critical 

loads  resulting  from  boundary  conditions   or  ratios  of  N   to  N  not 

o     i 

considered  in  tlie  present  work  maj-  satisfy  equation  (116)  within  tlie 
needed  radius  of  convergence.   It  would  t}icn  become  necessary  to  eithei' 
expand  the  solution  of  equation  (93)  about  some  other  point,  or  employ 
the  techniques  associated  with  analj'tic  continuatiori.   Tne  difficulties 
encountered  in  tlie  latter  appi'oacli  would  be  enormous. 


CHAPTER  VI 


COxXCLUSION 


The  px-esent  work  investigates  the  buckling  of  annular  sandwich 
panels.   Equilibrium  equations  and  boundary  conditions  satisfying 
continuity  requirements  were  derived  in  cartesian  coordinates,  using 
the  theorem  of  minimum  potential  energy.   These  equations  were  then 
transformed  into  polar  coordinates  througli  the  application  of  tensor 
analj'sis. 

Axisymmetric  buckling  being  assumed,  and  the  bending  rigidity 
of  the  faces  being  neglected,  the  equilibrium  equations  were  uncoupled 
by  using  a  n.odified  technique.   The  governing  equations  were  then 
compared  witli  existing  theories  for  single-layer  annular  panels  [11,12] 
and  circula]'  sandwich  panels  [10]. 

For  the  general  problem  of  an  annular  sandwich  panel  subjected 
to  unequal  inner  and  outer  compressive  loads,  and  constrained  by 
boundary  conditions  sim.ilar  to  those  emploj^ed  by  Olsson  [12],  a  power 
series  solution  was  obtained.   This  series  was  sliov.ii  to  possess  a 
radixis  of  con\'ergence  of  sufficient  magnitude.   Successive  appi'oxima- 
tions  were  then  computed,  and  a  graphical  solution  was  employed  for 
various  i-atios  of  outer  to  inner  compressive  loads.   Results  from  the 
f  if  til  approximation,  which  were  shown  to  be  within  2.2  per  cent  of  the 
exact  solution,  were  compared  with  those  obtained  from  earlier  theories 
[10,11,1?]. 
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The  present  work  represents  the  first  attempt  to  analyze  the 
stability  of  annular  sandwich  panels.   Further  extensions  of  the 
present  theory  raay  be  carried  out  bj'  including  the  effects  of  the 
bending  rigidity  of  the  faces  or  considering  boundarj'  conditions 
other  than  tliose  employed  here.   Furthermore,  continued  efforts  should 
be  directed  toward  obtaining  a  solution  to  the  unsymmetric  buckling 
pi'oblom.   In  tliis  way,  tlie  assumption  of  axisynunetric  buckling  could 
be  justified,  and  the  problems  associated  with  angular  dependent 
loading  functions  could  be  analj-zed. 


BIBLIOGRAPHY 


1.  F.  J.  Plant ema,  "Sandwich  Construction,"  Wiley  and  Sons,  New  York, 
1966. 

2.  S.  Timoshenko  and  S.  Woinowsky-Krieger ,  "Theory  of  Plates  and 
Shells,"  2nd  ed. ,  McGraw-Hill,  New  York,  1959. 

3.  S.  Timoshenko  and  J.  Gere,  "Theory  of  Elastic  Stability," 
2nd  ed.  ,  McGraw-Hill,  New  York,  1961. 

4.  N.  J,  Hoff,  "Bending  and  Buckling  of  Rectangular  Sandwich  Plates," 
NACA  TN  2225,  1950. 

5.  A.  C.  Eringen,  "Bending  and  Buckling  of  Rectangular  Sandwich 
Plates,"  Proc,  First  U.S.  Natl.  Cong.  Appl .  Mech. ,  1951,  pp.  381- 
390. 

6.  C.  C.  Chang  and  I.  K.  Ebcioglu,  "Elastic  Instability  of  Rectangu- 
lar Sandwich  Panel  of  Orthotropic  Core  with  Different  Face  Thick- 
nesses and  Materials,"  Transactions  of  the  American  Society  of 
Mechanical  Engineers,  J.  App.  Mech.,  Vol.  27,  No.  3,  September 
1960,  pp.  474-480. 

7.  S.  J.  Kim,  "Symmetric  and  Antisynmietric  Buckling  of  Sandwicii 
Panels,"  Doctoral  Dissertation,  Dept.  Eng.  Sci.  and  Mech., 
University  of  Florida,  1969. 

8.  E.  Reissner,  "Small  Bending  and  Stretching  of  Sandwicla-Type 
Shells,"  NACA  TN  1832,  1949, 

9.  M.  Zaid,  "Symmetrical  Bending  of  Circular  Sandwich  Plates,"  Proc, 
of  the  2nd  U.S.  Natl.  Cong.  Appl.  Mech.,  1951,  pp.  413-'422, 

10.  J.  C.  Huang  and  I.  K.  Ebcioglu,  "Circular  Sandwich  Plate  under 
Radial  Compression  and  Thermal  Gi'adient,"  AIAA  Journal,  Vol.  3, 
No.  6,  June   1965,  pp.  1146-1148. 

11.  E.  Meissner,  "Uber  das  Knicken  kreisringf ormiger  Scheiben," 
Schweiz.  Bauztg. ,  Vol.  101,  1933,  pp.  87-89. 

12.  G.  Olsson,  "Knicknng  der  Kreisringplatte  vcn  quadratisch 
veranderlicher  Steif igkeit ,"  Ingr.-Arch.,  Vol.  9,  1938,  pp.  205- 
214, 

73 


74 


13.  C.  C.  Chang  and  I.  K.  Ebcioglu,  "Thermoelastic  Behavior  of 
a  Simply  Supported  Sandwich  Panel  Under  Large  Temperature 
Gradient  and  Edge  Compression,"  J.  Aero.  Sci.,  Vol.  28,  No.  6, 
June  1961,  pp.  480-492. 

14.  I.  S.  Sokolnikoff,  "Mathematical  Theory  of  Elasticity,"  2nd  ed. , 
McGrav,--}Iill  ,  New  York,  1956. 

15.  Y.  C.  Fung,  "Foundations  of  Solid  Mechanics,"  Prentice-IIall , 
Englewood  Cliffs,  New  Jersey,  1965. 

16.  A.  J,  McConnell,  "Applications  of  Tensor  Analysis,"  Dover 
Publications,  Inc.,  New  York,  1957. 

17.  S.  Tim.oshenko,  "Strength  of  Materials,"  Part  II,  D.  van  Nostrand, 
Princeton,  New  Jersej',  1956. 

18.  F.  B.  Hildebranri,  "Advanced  Calculus  for  Applications,"  Prentice- 
Hall,  Englewood  Cliffs,  New  Jersey,  1962. 

19.  L.  R.  Ford,  "Differential  Equations,"  2nd  ed. ,  McGraw-Hill, 
New  York,  1955. 


BIOGRAPHICAL  SKETCH 

Amelio  John  ^mato  was  born  in  Newark,  New  Jersey,  on 
January  20,  1944.   He  was  graduated  from  Seton  Hall  Preparatory 
School  in  June,  1962.   In  June,  1966,  he  received  the  degree  of 
Bachelor  of  Science  in  Mechanical  Engineering  from  Newark  College 
of  Engineering  (New  Jersey). 

In  September  of  the  same  year  he  entered  tlie  Department 
of  Engineering  Science  and  Mechanics  at  the  Universitj'  of  Florida 
as  a  National  Defense  and  Education  Act,  Title  IV  Fellow.   Here, 
in  August,  1967,  he  received  tlie  degree  of  Master  of  Science  in 
Engineer jng  prior  to  pursuing  the  degree  of  Doctor  of  Philosophy. 


75 


This  dissertation  was  prepared  under  the  direction  of  the 
chairman  of  the  candidate' s  supervisory  committee  and  has  been 
approved  by  all  members  of  that  committee.   It  was  submitted  to  the 
Dean  of  the  College  of  Engineering  and  to  the  Graduate  Council,  and 
v/as  approved  as  partial  fulfillment  of  the  requirements  for  the  degree 
of  Doctor  of  Philosopliy. 


June,  1970 


Dean,  College  of  Engineering 


Supervisox-y  Committee: 


Dean,  Graduate  School 


>l7r^^'V  A    ^"^^/^.yk?^^ 


Aa.  ^A-  \Cvvx-»-w>e  \f 


tJLlMSd.^ 


